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ABSTRACT 


The  Surface  Structure,  Scattering  Losses  and 
Schottky  Barrier  Model  of  III-V  Compound  Semiconductors 

Professor  Bansang  W.  Lee 

The  surface  atomic  structure  of  GaP  (110)  face  has  been  deter¬ 
mined  by  low  energy  electron  diffraction  (LEED)  analysis.  LEED 
intensities  have  been  recorded  for  the  (110)  and  (111)  faces  of  GaP. 
An  analysis  of  the  structure  of  the  (110)  surface  of  GaP  is  performed 
by  comparing  LEED  Intensities  calculated  using  dynamical  multiple 
scattering  programs  with  those  measured  at  room  temperature.  The 
Intralayer  multiple  scattering  is  treated  exactly  using  direct  summa¬ 
tion  method,  while  for  the  interlayer  multiple  scattering,  the 
renormalized-forward-scattering  method  is  used.  The  resulting  best- 
fit  structure  suggests  that  the  P  atoms  protrude  from  the  GaP  (110) 
surface  plane,  whereas  the  Ga  atoms  are  rotated  inward  from  the  plane 
by  angle  of  27° ,  and  the  surface  layer  is  compressed  by  5Z.  The 
recorded  LEED  intensity  profiles  for  GaP  (111)  have  been  analyzed. 

The  result  suggests  that  a  good  Schottky  contact  would  be  relatively 
easy  to  obtain  on  the  GaP  (111)  face. 


Schottky  barriers  have  been  fabricated  on  n-type  GaP  (111) 
crystals.  Capacitance-voltage  and  current -voltage  measurements  were 
performed  to  analyze  the  electrical  properties  of  the  Schottky  diodes 


Auger  electron  spectroscopy  data  were  used  to  characterize  the  inter¬ 
face  chemical  composition.  A  metal-amorphous  film-semiconductor  (MAS) 
configuration  is  suggested  to  be  the  physical  structural  model  of  a 
real  Schottky  barrier.  Based  on  the  MAS  structural  model,  the  origin 
of  Fermi-level  pinning,  the  deviation  of  the  diode  factor  from  its 
ideal  value  and  the  transition  of  interface  electronic  behavior  be¬ 
tween  covalent  and  ionic  semiconductors  have  been  explained  and 
discussed . 

The  surface-scattering  losses  of  GaAs  waveguide  have  been  inves¬ 
tigated.  The  simple  theory  of  surface  scattering  and  the  effective 
index  approximation  were  combined  to  calculate  surface-scattering 
losses  of  a  rib  waveguide.  Different  degrees  of  surface  roughness 
were  prepared  by  polishing  the  GaAs  substrates  with  alumina  powders 
of  different  particle  sizes.  Laser  scattering  technique  was  then 
conducted  to  characterize  the  surface  geometry  of  these  substrates. 

The  result  shows  that  the  degrees  of  surface  roughness  of  the  sub¬ 
strates  were  proportional  to  the  grit  size  used  to  polish  them. 
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CHAPTER  1 


INTRODUCTION 


III-V  compound  semiconductors  are  used  widely  in  3olid  state 
electronic  devices  such  as  light  emitting  diodes  (LED) ,  Gunn  effect 
oscillators,  and  high  frequency  field  effect  transistors  on  account 
of  their  direct  and  relatively  large  bandgaps.  The  free  surface 
properties  of  III-V  compound  semiconductors  differ  profoundly  from 
those  of  the  elemental  semiconductors.  These  differences  Include 
surface  crystallography,  electronics  surface  states  and  surface 
chemical  activity.  A  special  role  has  been  played  by  GaAs  and  other 
III-V  compound  semiconductors,  which  appear  to  be  promising  in 
promoting  monolithic  integrated  optical  circuits.  All  optical 
function  must  be  obtained  in  a  single  substrate  in  monolithic 
integrated  optical  circuits.  The  optical  components  needed  to 
perform  complex  optical  functions  can  be  fabricated  in  a  single 
crystal  of  GaAs.  Most  of  the  III-V  compound  semiconductors  have 
large  electro-optic  effect  and  good  electrical  properties  for 
lasering,  wavegulding,  modulating,  switching  and  detecting  of  opti¬ 
cal  signals.  The  most  promising  substrate  material  known  to  date  is 
GaAs.  However,  the  optical  transmission  losses  in  GaAs  waveguides 
cannot  be  reduced  to  an  acceptable  level  by  the  present  technologies 
for  material  preparation.  A  detailed  understanding  of  the  magnitude 
of  the  propagation  loss  in  semiconductor  waveguides  is  needed.  Rec¬ 
ognizing  the  useful  application  and  the  urgent  needs  in  the  field 
of  optical  electronics,  it  is  necessary  to  gain  a  quantitative 


understanding  of  the  free  surface  and  interface  properties  of  the 
III-V  compound  semiconductors  in  order  to  develop  a  fabrication 
technology  of  the  same  quality  that  exists  for  silicon-based  devices. 

One  of  the  most  fundamental  types  of  information  needed  about 
solid  surfaces  concerns  the  relative  atomic  positions.  The  geomet¬ 
rical  arrangement  of  surface  atoms  influences  most  physical  and 
chemical  properties  of  surfaces,  such  as  electronic  surface  states, 
oxidation,  contact  potentials  and  crystal  growth.  It  has  been  found 
that  low  energy  electron  diffraction  (LEED)  is  a  successful  tool  to 

determine  the  relative  atomic  positions  at  surfaces.  Incident  elec- 

o 

crons  in  the  energy  range  10-400ev  penetrate  not  more  than  5-iOA 
into  the  surface  region.  The  wave  interference  basic  to  LEED  con¬ 
tains  an  abundance  of  information  about  the  geometrical  configuration 
of  atoms  at  surfaces.  LEED  can  be  used  quantitatively  to  determine 
the  atomic  structure  of  single  crystal. 

In  this  thesis,  LEED  and  Auger  electron  spectroscopy  are  used 
to  investigate  the  free  surface  properties  of  GaP.  Modern  surface 
analytical  techniques  and  standard  I-V  (current-voltage)  and  C-V 
(capacitance-voltage)  techniques  are  applied  to  characterize  the 
physical  mechanism  for  the  metal-semiconductor  interaction  and  for¬ 
mation  of  the  Schottky  barriers.  In  regard  to  the  waveguide 
transmission  losses  and  surface  roughness,  a  quantitative  character¬ 
ization  of  the  surface  "smoothness"  of  GaAs  substrate  has  been  made 
before  the  substrate  is  fabricated  into  a  waveguide,  and  this  quali¬ 
tative  information  of  the  surface  will  be  used  to  develop  the 


correlation  between  the  surface  roughness  and  scattering  losses  of 
the  waveguide. 

The  theoretical  formulism  for  LEED  intensity  analysis  is  re¬ 
viewed  in  the  first  part  of  the  next  chapter.  The  notation  of 
surface  structure  and  the  physical  origin  of  the  LEED  patterns  are 
first  described  and  then  the  physical  model  for  a  one-electron 
schrodinger  equation  is  constructed  from  which  the  multiple  scat¬ 
tering  theory  referred  to  as  dynamical  diffraction  theory  can  be 
developed.  The  dynamical  formulation  used  in  our  theoretical  cal¬ 
culation  is  introduced  in  the  last  section.  The  effects  of  surface 
atomic  thermal  vibration  (electron-phonon  interaction)  and  inelastic 
damping  which  removes  electrons  from  the  flux  of  elastically  scat¬ 
tered  electrons  in  the  scattering  process,  are  also  included  in  the 
theoretical  calculation  to  Improve  the  agreement  of  the  calculated 
LEED  intensities. 

The  theoretical  calculation  of  surface  scattering  theory  is 
reviewed  in  the  second  part  of  the  next  chapter.  The  surface-scat¬ 
tering  of  a  symmetric  slab  waveguide  calculated  on  the  basis  of 
radiation  mode  is  first  described  and  then  a  simple  theory  based  on 
the  Rayleigh  criterion  is  introduced  to  calculate  the  surface  scat¬ 
tering  losses  for  an  asymmetric  slab  waveguide.  Since  it  is 
impossible  to  measure  the  scattering  losses  of  the  two  asymmetric 
surfaces  separately,  all  the  surface  properties  were  lumped  into  a 
single  parameter  in  the  theory.  In  integrated  optics,  however, 
rectangular  strips  of  dielectric  material  embedded  in  other  dielec- 


cries  are  used.  The  theoretical  formulation  of  surface-scattering 
losses  for  a  rectangular  waveguide  is  derived  in  the  last  section. 

An  exact  analytical  treatment  of  this  problem  is  impossible,  there¬ 
fore,  an  approximated  analytical  approach  is  used  together  with  the 
simple  theory  in  solving  it. 

The  experimental  methods  are  described  in  Chapter  Three.  The 
physical  origin  of  Auger  electron  spectroscopy  (AES)  and  the  exper¬ 
imental  set  up  for  obtaining  the  Auger  spectra  are  described.  The 
technique  used  to  characterize  the  surface  geometry  of  the  waveguide 
is  introduced.  The  method  used  in  forming  a  rib  waveguide  for  in¬ 
vestigating  the  surface-scattering  losses  of  III-V  compound 
semiconductor  waveguides  is  suggested.  Finally,  all  the  electrical 
measurements  with  actual  data  handling  are  described. 

In  Chapter  Four,  the  dynamical  programs  formulated  in  Chapter 
Two  are  applied  to  the  GaF  (110)  surface  for  surface  atomic  structure 
determination.  LEED  intensities  recorded  experimentally  with  normal 
incident  primary  electron  beams  at  both  room  temperature  and  low 
temperature  were  used  in  analyzing  the  surface  structure.  The  atomic 
structure  information  such  as  layer  spacing  and  atomic  positions  are 
quantitatively  obtained  through  the  analyses  of  the  LEED  patterns 
and  the  intensity  profiles.  The  LEED  patterns  and  experimental 
intensities  for  the  GaP  (111)  face  are  analyzed,  and  the  surface 
structure  of  GaP  (110)  face  is  compared  with  that  of  GaAs  (110) 


In  Chapter  Five,  the  physical  mechanism  for  the  metal-semicon¬ 
ductor  interaction  and  forming  of  the  Schottky  barriers  are  discussed. 
As  suggested  in  Chapter  Four,  the  GaP  (Ill)  face  was  chosen  to  be  the 
semiconductor  surface  for  Schottky  barriers  studies.  A  macroscopic 
structural  model  of  the  real  Schottky  barrier  is  suggested.  On  the 
basis  of  the  structural  model,  an  electrical  equivalent  model  is 
derived  to  address  the  fundamental  behaviors  of  a  real  Schottky 
diode,  including  the  origin  of  Fermi-level  pinning,  the  barrier 
heights,  and  the  deviation  of  the  diode  factor.  Schottky  barriers 
were  fabricated  on  GaP  (111)  face  to  demonstrate  the  structural 
model.  Electrical  measurements  described  in  Chapter  Three  were 
performed  to  analyze  the  electrical  properties  of  the  diodes.  AES 
were  used  to  correlate  the  chemical  composition  of  the  metal- 
semiconductor  Interface  to  the  structure. 

The  surface-scattering  losses  of  GaAs  waveguide  are  investi¬ 
gated  in  Chapter  Six.  The  simple  theory  of  surface  scattering 
introduced  in  Chapter  Two,  together  with  the  effective  index  approx¬ 
imation,  is  used  to  calculate  the  surface-scattering  losses  of  a  rib 
waveguide.  The  experimental  methods  to  prepare  different  degrees 
of  surface  "smoothness"  of  GaAs  substrates  are  described.  The  re¬ 
sults  of  quantitative  characterization  of  the  surface  geometry  of 
the  substrates  are  reported.  The  surface  geometric  information 
combined  with  the  signal  losses  measured  after  the  substrate  is 
fabricated  into  a  rib  waveguide  will  be  used  to  develop  the  corre¬ 
lation  between  the  surface  "smoothness"  and  the  surface  scattering 
parameters . 


CHAPTER  2 


TECHNICAL  BACKGROUND 

2.1  Review  of  Theoretical  Formulation  for  Determination  of  the  Atomic 
Structure  of  Solid  Surface  by  Low  Energy  Electron  Intensity  Analysis 

2.1.1  Introduction 

It  has  been  known  for  a  very  long  time  that  low  energy  electron 
diffraction  (LEED)  is  a  powerful  tool  in  the  study  of  the  atomic 
structure  of  solid  surfaces.  However,  It  is  only  since  the  achieve¬ 
ment  of  ultrahlgh  vacuum  (UHV)  technology  that  LEED  has  been  widely 
used  and  studied. 

By  measuring  and  analyzing  coherent  elastically  backscattered 
electrons  from  a  crystal  surface,  LEED  can  provide  the  geometric 
structure  of  a  free-surface  (solid-vacuum  interface)  of  a  solid, 
namely,  the  unit  mesh  symmetry,  the  unit  mesh  constants,  the  atomic 
position  of  the  same  or  different  kinds  of  atoms  in  the  unit  mesh, 
and  furthermore,  the  spatial  distance  between  the  surface  atomic 
layers,  the  shift  distance  of  unit  mesh  relative  to  their  substrates, 
and  the  relationship  of  the  observed  surface  structure  to  the  nature 
of  chemical  bonds  of  the  substrate  surface. 
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2.1.2  Notation  for  Surface  Structure 

Surfaces  are  usually  named  after  the  planes  of  bulk  crystal  to 
which  they  are  parallel,  e.g.,  (100),  (110),  etc.  Their  two-dimen¬ 
sional  structure  can  be  characterized  in  terms  of  a  unit  mesh,  the 
smallest  unit  that  can  be  repeated  indefinitely  to  build  up  the 
whole  surface.  The  two  sides  of  the  unit  mesh  will  be  designated  by 
vectors  a  and  b  as  shown  in  Fig.  2.1.  The  origin  of  any  mesh  can  be 
given 


relative  to  the  origin  of  a  reference  unit  cell.  1  and  m  are  inte¬ 
gers.  In  the  two-dimensional  case,  there  are  five  kinds  of  Bravaise 
lattice  as  shown  in  Fig.  2.2(a).  The  corresponding  reciprocal  lat¬ 
tice  are  shown  in  Fig.  2.2(b).  The  unit  mesh  vectors  a*  and  b*  of 
the  reciprocal  lattice  are  given  by^ 

a*  -  2*(a  b  -b  a  )_I(b  ,-b  ) - (2.2) 

x  y  x  y  y  x 

b*  «  2*(a  b  -b  a  )"l(-a  ,a  ) - (2.3) 

x  y  x  y  y  x 

where  a^  and  b^  are  the  x  components  of  a  and  b,  respectively,  and 
a  and  b  are  the  y  components  of  a  and  b,  respectively. 


A 
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If  a  reconstructed  surface  is  to  be  described  in  terms  of  the 

bulk  exposed  plane  or  if  a  monolayer  of  some  surface  deposit  has  a 

different  unit  mesh  from  that  of  the  solid  surface  upon  which  it  is 

2 

adsorbed,  a  simple  notation  devised  by  Wood  is  used.  Two  examples 
of  this  notation  are  illustrated  in  Fig.  2.3.  The  unit  meshes  of 
the  surface  and  bulk  are  related  by  an  expression  of  the  form 
(ag/a^  x  bg/b^)R.  The  subscripts  s  and  b  are  used  for  a  surface 
mesh  and  for  a  bulk  exposed  plane,  respectively.  If  the  deposit 
mesh  is  not  rotated  with  respect  to  the  substrate  (or  bulk)  mesh, 
then  the  rotation  R  is  simply  dropped  from  the  notation.  If  the 
deposit  mesh  is  a  centered  mesh,  then  a  lower-case  c  will  be  used 
before  the  expression. 


[n] 


M^.NiU10)  c  (2X2) 

substrate  substrate  centered  .  surface/  deposit 
surface  mesh  substrate 
nor  ma  I  mesh  ratio 

(a) 


•  O  Cfr-fcft  •  o  o  o 


O'  L 

©•--■©  o  •  o  o 


0  0  °  • 

O  ,0  o  >  o  o  o 

/  / 

<  O  0X  O  #00 


os  V  o  o  o 


o  o  o 
•  o  o 


e.g.  Pt(lOO)  ( /2X/T)R45°-0 


Fig.  2.3  Two  examples  of  Wood's  notation  for  surface  structures; 

(a)  Ni  (110)  face  exposed  to  oxygen,  (b)  P  (100)  face 
exposed  to  oxygen. 


2.1.3 


The  Diffraction  Pattern 


Assume  a  mono-energetic,  collimated  beam  of  electrons  incident 
normally  on  a  perfect  two-dimensional  lattice.  Because  of  the  peri¬ 
odicity  along  the  direction  of  surface  normal  is  lost  in  a  two-dimen¬ 
sional  case,  the  constructed  interference  of  scattered  waves  cannot 
occur  along  this  direction.  This  condition  leads  to  the  possibility 
of  diffracted  beam  occurring  at  all  incident  energies  and  hence  a 
diffraction  pattern  can  be  observed  at  all  energies. 

The  unit  mesh  vectors  a*  and  b*  of  the  reciprocal  lattice  are 
related  to  the  unit  mesh  vectors  a  and  b  of  the  real  lattice  via 
Eqs.  (2.2)  and  (2.3).  The  energy  E  of  the  incident  electron  beams 
is  given  by 

E  -  t(2k2  /2m - (2.4) 

with  k  -  2ir/X  and  fi  -  h/(2ir) - (2.5) 

where  h  is  Planck's  constant,  m  is  the  mass  of  electron  and  X  is  the 
wavelength. 

The  Ewald  sphere  construction  can  then  be  applied  to  the  two- 
dimensional  diffraction  problem  as  shown  in  Fig.  2.4.  Because  of 
the  loss  of  periodicity  in  one  dimension,  the  reciprocal-lattice 
rods  can  be  labelled  with  two  Miller  indices  h  and  k,  and  a  general 
reciprocal  lattice  vector  g^  lies  in  the  plane  of  the  surface  and 


is  given  by 


ghk-ha*+  kb* - —(2.6) 

Because  every  point  on  the  reciprocal-lattice  rods  is  a  reciprocal 
lattice  point,  diffraction  occurs  at  the  points  where  Evald  sphere 
cuts  the  reciprocal  rods  and  the  diffracted  beam  can  be  labelled 
with  the  Miller  indices  (hk)  of  the  rod  causing  it.  The  wave  vector 
of  the  diffracted  beam  can  be  characterized  by 

- (2-7> 

Suppose  the  backscattered  electrons  are  detected  on  a  fluorescent 
screen.  These  points  are  visible  on  the  screen  as  bright  spots 
which  constitute  the  electron  diffraction  pattern  of  the  two-dimen¬ 
sional  lattice.  As  shown  in  Fig.  2.4,  the  pattern  provides  us  with 
a  direct  picture  of  the  reciprocal  lattice  and  an  image  of  crystal 
periodicity.  Hence,  the  diffraction  pattern  immediately  reveals 
the  unit  mesh  symmetry  and  dimension  ratio. 

As  the  energy  of  the  incident  electron  changes,  the  magnitude 
of  the  radius  of  the  Evald  sphere  changes.  So  long  as  the  radius 
is  large  enough  to  let  the  sphere  intersect  the  reciprocal-lattice 
rods,  the  spots  of  the  diffraction  always  exist  with  constant  in¬ 


tensity. 
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In  nature,  the  two-dimensional  case  does  not  exist  because  the 

o 

incident  electrons  penetrate  a  distance  of  8  to  12  A  into  the  bulk 
of  the  crystal,  they  go  sufficiently  deep  to  be  sensible  of  the 
periodicity  along  the  direction  of  the  surface  normal.  Constructed 
interference  of  scattered  waves  occurs  in  this  direction.  Further 
information  about  the  deposition  of  contents  of  the  cell,  such  as 
separation  between  planes  of  atoms  or  location  of  absorbers  relative 
to  the  substrate,  which  cannot  be  provided  by  the  diffraction  pattern, 
is  needed.  This  information  is  contained  in  intensities  of  the  back- 
scattered  beams  measured  as  functions  of  incident  energy  or  angle  of 
incidence.  Intensity/ energy  curves  are  the  usual  form  for  presenting 
data.  The  theoretical  formulation  of  LEED  intensity  curves  will  be 
discussed  in  section  2.5. 


< 


2.4  Generation  of  LEED  pattern  from^a  two-dimension  lattice, 
electron  beam  with  wave  vector  k0  is  normally  incident 
upon  the  crystal. 


2.1.4  Physical  Model  Used  in  LEED  Intensity  Calculation 

The  physical  model  of  a  crystal  surface  is  assumed  to  be  a 
defect-free,  single  crystal  with  an  atomically  flat,  absorbent-free 
surface  of  at  least  several  square  millimeters  in  area.  All  theo¬ 
retical  formulations  and  analyses  of  surface  structure  in  LEED  are 
based  on  this  assumption.  Developments  in  LEED  over  the  past 
several  years  have  demonstrated  that  this  highly  idealized  surface 
gives  valuable  guidance  in  the  interpretation  of  experimental  results 

The  first  step  of  LEED  intensity  calculation  and  analysis  is  to 
predict  the  peak  positions  in  the  intensity  vs.  energy  curves  of  a 
certain  spot  in  the  diffraction  pattern.  Bragg  Theory  was  combined 
with  electron  scattering  theory  to  develop  a  formulation  of  LEED 
intensity  some  twenty  years  ago.  This  method  is  called  the  "kine- 
matical  model".  The  assumption  of  this  formulation  is  that  the 
Incident  electrons  are  scattered  elastically  once  and  only  once  from 
the  atoms  of  the  crystal.  This  kinematic  computation  of  LEED  inten¬ 
sity  only  gives  some  quantitative  understanding  of  the  crystal  sur¬ 
face.  It  can  only  predict  the  positions  of  Bragg  peaks.  Actually, 
in  an  experimental  intensity  curve,  there  is  much  more  structure 
than  those  predicted  by  the  kinematic  theory.  "Secondary"  peaks, 
"shoulders"  and  other  peak  shape  anomalies  in  measured  intensity 
curves  can  only  be  formulated  by  multiple  scattering  of  the  incident 
electrons  from  the  crystal  lattices.  Multiple  scattering  events 
thus  must  be  used  in  the  formulation  of  LEED  theory  to  obtain  a  good 
description  of  the  intensity  curve  and  to  gain  a  quantitative 


understanding  of  the  crystal  surface.  The  multiple  scattering  for¬ 
mulation  is  referred  as  the  "dynamical  model"  of  LEED  theory. 
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There  have  already  been  abundant  articles  in  the  literature  ’  * 
reviewing  the  kinematic  approach  to  the  theory  of  LEED  in  detail,  so 

it  will  not  be  discussed.  The  dynamical  formulation  constructed  from 

a  one-electron  Schrodinger  equation  will  be  developed. 


As  the  incident  electrons  approach  the  crystal,  the  most  impor¬ 
tant  interaction  is  the  elastic  scattering  of  the  incident  electron 
by  ion-cores  comprising  the  lattice.  Let  us  consider  a  system  con¬ 
sisting  of  a  solid  rigid  lattice  and  a  probing  electron.  We  should 
be  solving  the  complete  Schrodinger  equation  that  includes  all  elec¬ 
trons  in  the  total  wavefunction  <5. 


+ 


N 


(2.8) 


Here,  the  electrons,  labelled  by  the  index  j,  have  positions  r^  and 
spin  coordinates  s^ ,  the  nuclei  are  labelled  by  a,  and  have  atomic 

number  Z  and  positions  R  ,  an  incident  electron  has  position  ro 

ci  a 

and  spin  coordinate  s<j.  E£  is  the  total  energy  which  includes  the 

core  levels  as  well  as  that  of  the  incident  electron.  V  is  the 

s 

potential  due  to  screening  charge.  By  Pauli  principle,  the  wave- 
function  <t  is  antisymmetric  under  exchange  of  particles  because  elec¬ 
trons  are  Fermions.  We  can  write  down  <J>  as  a  product  of  core-state 


wavef unctions . 


$(r0,30,  — ,rN,sN)» 


E  (-1)  (Kro  ,so)^i  (r  ,s  )■ 
P  J  J 


-^(rk,sk)-(2.9) 


<j>  is  the  wavef  unction  of  incident  electron.  The  sum  is  over  permuta¬ 
tions  of  electron  coordinates.  There  are  OW-1)!  possible  permutation 
Since  the  core-state  part  of  4  is  know  already,  it  should  be  elimi¬ 
nated  from  Eq.  (2.8)  to  simplify  the  problem. 


Assume  the  individual  wavef unctions  are  orthonormal 


E  JV*  (r,s)  H'j  (r,s)d3  r  »  6^ 


E  Pt*  (r,s)  <fi  (r,s)d3  r  -  0 


■(2.10) 


and  define  a  function 


X ( r :x ,  s L , - rN ,  sN)  -  f  L  ( r x  ,  s x )  ¥ 2  ( r 2 ,  s2) - ^ ,  s^) - (2.11) 


★ 

Multiply  Eq.  (2.8)  by  X  ,  and  then  integrate  it  over  space  and  sum  it 
over  spin  coordinates.  The  only  term  remaining  on  the  right-hand 
side  of  the  equation  is  that  in  which  the  coordinates  have  not  been 
permuted,  leaving  just  Et4(ro,s0)  after  integration  and  summation. 

The  left-hand  side  of  the  equation  can  be  split  into  three  parts 


i  diJ1— d3?/  c  d-nW3-^/ 1^4  |  +  v.dj) 

S,  „  1-1  "  j  ., 


N  N 


+  Z  Z  l/lr.-r.l)  +  (  -«2/2m7?  -  ZZ  /|S .  -r0  |  +  V  (r0>) 
1-1  i-j+1  J  1  ot  a  a  s 


j-1  i-j+1 


+  (Zl/| ro-r. | )  )  $ 


The  first  term  is  zero  unless  the  permutation  in  <p  are  those  leaving 
<*>-<J>(ro  ,so) .  From  the  first  term  of  the  above  equation  we  obtain 
EcC^(ro,So),  where  E  is  the  total  energy  associated  with  the  core¬ 
states.  The  second  term  has  non-zero  value  only  if  the  coordinates 
of  Ohave  not  been  permuted,  i.e., 


(  -ti/2m  -ZZa/|Ra-^oi  +  V3(r0)  )  ♦(r0,s0) 


Finally,  the  third  term  has  non-zero  contribution  only  from  the  first 
component  of  4>aad  from  those  permutations  that  exchange  coordinates 
of  <Kro,s0) 


N 

Z  <Kro,so)/Z'i,*<rj,Sj)  ^  (r^  ,3^)  /  j  r0-r^  |  d3^  - 
j— 1  sj 


Z  f.  (ro  .So)/  Z'Mr  ,s  )  <t>(r  ,s  )/|ra-r .  Id3^ 

j»l  J  sj  J  J  J  JJ  JJ 


The  complete  equation  is  reduced  to 


(-«2/2mV?-  ZZ  /|R  -raj  +  V  (r„>  +  Zf  I|V. (r. ,s.) | 2/ 

a  a  a  s  j-i  3j  J  3  J 

N 

|r0-r [d3?  .)<j>(r0  ,s0)-  Z  (/  E  4f*(r  ,s .)4>(ra ,30) / 1 r0-r  | cT^r  ) 

J  J  gj  j  j  J  j  J 

¥  (r^.s  )  »(Et-Ect)  4>(r9t3o)-  E<J>(r0fs0) - (2.12) 


The  strong  resemblance  of  Eq.  (2.12)  to  a  simple  one-particle 
Schrodlnger  equation  enables  us  to  connect  the  usual  one-particle 
theories  of  diffraction  with  many-body  problems.  In  essence,  the 
usual  theories  of  diffraction  are  based  on  this  equation.  The  second 
term  in  Eq.  (2.12)  arises  because  of  antisymmetry  under  exchanges  of 
electrons,  for  no  two  electrons  can  appear  at  the  same  place  at  the 
same  time. 


Besides  the  electron  ion-core  interaction  there  are  other  con¬ 
tributions  during  the  scattering  of  the  electron  from  the  crystal. 
They  are  from  the  conduction  or  valence  electrons,  and  at  finite 
temperatures  from  phonons  as  well.  Electron  and  electron  inter¬ 
actions  are  strong  in  the  energy  range  of  interest  (lOev  and  300ev) . 
They  excite  surface  and  bulk  plasmons  and  are  responsible  for  energy 
losses  during  inelastic  scattering.  A  complete  detailed  microscopic 
formulation  and  analysis  of  this  interaction  is  beyond  the  existing 
capacity  of  physical  theory.  The  modification  is  generally  done  by 
adding  an  empirical  complex  "optical"  potential  Vq  in  the  Schrodin- 
ger  equation  (2.12)  for  inelastic  scattering. 


V  +  iV  ; 
or  oi 


•(2.13) 
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The  real  pare,  V  ,  of  Che  optical  potential  (also  known  as  the  inner 
or 

potential)  is  contributed  by  the  exchange  interaction  of  incident 
electron  with  the  conduction  or  valence  electrons  of  the  crystal. 

The  imaginary  part,  Vq^,  of  the  optical  potential  is  contributed  by 
the  inelastic  collision.  Both  V  and  are  treated  as  adjustable 
parameters  whose  values  are  selected  during  a  calculation  to  achieve 
the  best  agreement  with  experimental  intensity  profiles.  When  the 
temperature  of  a  crystal  is  raised,  two  changes  become  apparent  in 
the  diffraction  pattern:  firstly,  intensities  of  individual  diffrac¬ 
tion  spots  decrease  and  secondly,  the  intensities  of  the  background 
collected  between  spots  grow  stronger.  It  is  the  thermally  induced 
vibration  of  atoms  chat  are  responsible  for  these  phenomena.  The 
theoretical  treatment  of  the  temperature  effect  in  LEED  formation 
will  be  discussed  In  the  next  section. 


2.1.5 


Dynamical  Formulation  of  LEED 
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Several  methods  have  been  proposed  for  solving  Schrodinger 
equation  of  multiple  scattering.  One  of  the  methods*’  ^  mostly  used 
in  LEED,  is  to  formulate  the  dynamic  theory  from  the  multiple  scat¬ 
tering  point  of  view.  The  basic  viewpoint  of  multiple  scattering 
theory  is  that  the  interaction  between  the  incident  electron  beam 
and  scattering  centers  is  visualized  in  terms  of  a  succession  of 
scattering  effects.  Based  on  this  viewpoint  one  can  express  the 
entire  complex  scattering  problem  in  terms  of  the  individual  scat- 
terers.  A  multiple  scattered  electron  from  the  crystal  ion-core 
potential  can  be  viewed  as  composed  of  sum  over  paths  of  an  electron 
undergoing  individual  scattering  events  from  each  separated  ion-core 
potential.  The  whole  crystal  is  assumed  to  be  constructed  of  the 
largest  non-overlapping  spheres  of  each  atom.  The  potential  inside 
the  sphere  will  be  spherically  symmetric,  and  the  potential  outside 
will  be  constant.  This  is  called  the  "muffin  tin  model"  of  the 
crystal  potential.  The  radius  of  the  potential  sphere  is  called  the 
muffin-tin  radius  of  the  specific  atomic  species.  The  incident  elec¬ 
trons  are  assumed  to  be  scattered  elastically  by  the  individual 
potential  sphere  only  and  propagating  freely  in  the  spaces  between 
the  spheres.  Under  these  assumptions,  the  individual  scattering  of 
electron  from  the  sphere  can  be  characterized  by  a  phase  shift  of 
scattered  wave  relative  to  the  incident  wave.  To  implement  this 
procedure  the  incoming  plane  wave  is  expanded  in  spherical  harmonic 
partial  wave  components.  The  computation  of  I'th  partial  wave's 
phase  shift  <5^  (£-0,1,2--)  can  be  found  in  standard  text  of  quantum 


mechanics  .  Computer  programs  for  calculation  6^  are  available  in 
the  literature^.  Thus,  the  details  of  8^  computation  will  not  be 
discussed.  There  are  many  equivalent  mathematical  formulism  for 
solving  the  muliptle  scattering  equations,  and  there  are  correspon¬ 
dingly  many  computer  programs  to  Implement  these  various  schemes. 

The  computer  programs  used  in  our  dynamical  LEE D  intensity  calculation 
were  supplied  by  Tong,  et  al.  and  the  formulism  they  used  for  solving 
the  multiple  scattering  is  known  as  the  layer  method*,  in  which  the 
bulk  crystal  and  adsorbate  system  is  replaced  by  a  finite  slab  made 
by  layers  of  atoms  parallel  to  the  surface  of  the  crystal.  Multiple 
scattering  is  calculated  first  among  the  atoms  of  a  single  layer, 
and  subsequently  among  the  stack  of  layers  which  form  the  complete 
crystal.  The  method  will  be  described  as  follows: 

Assume  an  electron  beam  is  incident  on  the  crystal.  The  surface 
model  of  the  crystal  is  defined  in  terms  of  a  right-hand  coordinate 
axes,  two  of  which  (x,y)  lie  in  the  plane  of  the  layer,  parallel  to 
the  sufrace  with  the  third  (z)  pointing  into  the  crystal  and  defining 
the  positive  direction  for  wave  propagation.  Each  layer  of  a  stack 
consists  of  a  periodic  array  of  atoms  represented  for  scattering 
purposes  by  spherically  non-overlapping  muffin-tin  potentials.  The 
scattering  properties  of  each  atomic  potential  are  described  by  a  set 
of  phase  shifts,  and  the  constant  potential  outside  the  muffin-tin 
spheres  is  taken  as  the  zero  level  for  the  calculation.  In  general, 
the  position  vector  of  an  atom  within  a  layer  is  written  as 


r  is  the  location  of  the  a'th  atom  in  p'th  unit  cell. 

P  K 

r  -  la  +  mb - (2.15) 

P 

where  1  and  m  are  integers,  and  a  and  b  are  the  unit  mesh  vectors 
within  the  layer.  Eq.  (2.15)  restricts  the  number  of  atoms  in  the 
unit  cell  to  one  (ct*A)  or  two  (a*A,B) .  The  displacement  vector  s^ 
from  A-atom  to  B-atom,  may  have  a  non-zero  perpendicular  component, 
giving  a  non-coplanar  or  "rumpled"  layer.  The  incident  electron 
beam  is  given  in  terms  of  its  energy  E  in  the  vacuum,  and  the  two 
direction  angles,  e.g.,  the  polar  angle  9  measured  with  respect  to 
the  surface  normal  (+z)  ,  and  the  azimuthal  angle  <t>  measured  with 
respect  to  the  direction  chosen  as  (-he)  in  the  plane  of  the  surface. 
The  components  of  momentum  of  the  Incident  beam  parallel  to  the  sur¬ 
face  are  given  by 

ksin9cos0,  ksin0sin$ — - - (2.16) 

with  k-|k|  -(2mE/*V/2 - (2.17) 

If  a  plane  wave 

<t+  "  Aexp ( ik** r) - (2.18) 

is  incident  on  the  crystal,  the  2-d  periodicity  means  that  reflected 
waves  must  be  comprised  of  a  set  of  discrete  beams 


r  -  (  (kx  +  g^) ,  (ky  +  gy) ,  -(2mE/h2 

-2Vo  -  (V  gx)2  -  (ky  +  gy)V/2>-(2.20) 

where  g  is  a  reciprocal  vector  of  2-d  mesh,  Vq  is  the  optical  poten¬ 
tial  of  the  crystal,  and  E  is  the  energy  of  the  incident  plane  wave. 
We  choose  the  +z  direction  to  be  pointing  inward  into  the  crystal 
surface. 


Consider  a  set  of  plane  waves 


<P~  *  5  a-*  exp(±  *C*r> 
g  g  g 


■(2.21) 


incident  from  the  +z  and  -z  directions  on  a  layer.  It  is  more  con¬ 
venient  to  re-express  the  plane  waves  as  a  sum  of  spherical  partial 
waves.  The  spherical  wave  expansion  of  Eq.  (2.21)  about  the  a’th  atom 
in  the  p'th  unit  cell  becomes 


(T  -  £  a*  exp(i  k~  (r-r®))exp(i  k^  *r) 
g  g  g  P  g 


g£  Al£  Ylm(n(*'^))exp(iY*p) - (2>22) 


with 


Ao) 

*ma 


l  C  4iri*  (-1)“  T4  (Q(k*  ))exp(ik"  «s  ) 
g  g  g  g 


- (2.23) 


where  is  the  amplitude  of  (l,m)  angular  momentum  component  of 

the  incident  wave  at  the  a'th  ion-core  in  the  unit  mesh,  and 

fc-lk4!  =■  (2mE/h2-2V  )1/2 - (2.24) 

o 

s 


Q(r-rp)  stands  for  the  angular  coordinates  of  r-r^.  j  is  spherical 
Bessel  function.  is  the  spherical  harmonics  function  of  angular 

variation  of  the  wavefunction,  and  1  and  m  are  angular  quantum  numbers. 


For  multiple  scattering,  the  total  amplitude  of  the  waves  inci- 

ore 
(s) 


dent  on  the  a'th  ion-core  in  the  unit  cell,  Ap^,  will  include  a 


further  contribution  A^^  ,  representing  the  amplitude  scattered  from 
other  ion-cores,  which  depends  on  A..  Therefore,  these  quantities 
must  be  determined  by  self  consistent  iteration* . 


InCL 


A<«>  +  aJS) 
£ma  £ma 


■(2.25) 


(s) 

■rith  A^^  ,ima 


- (2.26) 


where  X.-  .  is  the  multiple  scattering  matrix  connecting  the 

Jt  t&  D  *  ~  m  Qt 

a'th  and  3'th  unit  meshes,  combining  Eqs.  (2.25)  and  (2.26)  gives 


-1 


A2ma  “  Z'm'SAl>m'S(  (I“X)  }  Z'm'8,lma“'(2*27) 


which  is  the  amplitude  of  the  scattered  beam. 


g  g  g  g  g 


where  M~  "  are  Che  scattering  matrices  of  the  layer.  The  sign  con- 

SS  -  + 

vention  is  such  that  M  describes  scattering  of  a  beam  g  incident 

g'Z  ^ 

in  positive  (+)  direction,  to  beam  %'  reflected  back  in  the  negative 
(-)  direction. 


<;  -  i-S-6 » 

g  g  g  z  g 


exp(ik;.sB))  ta(i\(!2(k;  » 

g  g 


exp(-ik~^*s^)  exp(id“  )sin6“) 
g'  * 


- (2.29) 


where  A  is  the  area  of  the  unit  cell,  and  6“  is  the  £th  wave  phase 
shift  of  the  atom  in  subset  a. 


This  completed  the  formulation  of  multiple  scattering  within  a 
single  layer.  The  multiple  scattering  subsequently  among  the  stack 
of  layers  will  be  examined  in  the  next  section. 


The  division  of  the  crystal  into  layers  is  illustrated  in 
Fig.  2.5,  which  shows  a  cross-sectio *  normal  to  the  surface  through 
the  atom  cores  in  successive  layers.  The  layers  are  labelled  i*0,l, 

2 -  from  the  vacuum  into  the  crystal,  i»l  denotes  the  surface  layer 

and  i»0  denotes  the  position  of  the  surface  potential  barrier,  layer 
displacement  d.  are  measured  into  the  crvstal  (+z-direction)  and 


Pig.  2.5  Cross  section  normal  to  the  surface,  showing  the  potential  inside  the  crystal 
and  the  arrangement  of  the  layers. 


measured  from  an  A-atom  site  in  the  i'th  layer  to  an  A-atom  site  in 
i+l'th  layer.  The  displacement  of  the  surface  potential  barrier  from 
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the  first  layer  of  atom  d^  defines  the  physical  position  of  the  sur¬ 
face.  It  is  a  adjustable  parameter.  The  layers  are  separated  by 
constant  potential  regions  and  in  these  regions  the  wavefield  can  be 
represented  by  a  superposition  of  plane  waves,  e.g., 

<j>i(r)»E  a+_^  exp(ik^*r)+a  _^exp(ik^*r) - (2.30) 

g  ig  g  ig  g 

where  r  is  referred  to  a  point  halfway  betweer  1ayers.  Scattering 
of  such  a  wavefield  through  the  i'th  layer  is  governed  by  the  equation, 


+  _+  ++  _+  +  +  +  —  —  - 

Ai  "  Pi  Ti  Pi-lAi-l  +  Pi  Ri  Pi  Ai 


Ai-i  -  pi-i  ptiAt-i + pi-iTr  pi\ 


■(2.31) 


where  a  vector  notation  for  the  beam  amplitudes  has  been  used,  e.g.. 


i 


A~  «  (a  'J) . 
1  ig 


The  matrices 


T^  =■  I  +  M** 
i  i 


I  +  M 

i 


< 


-(2.32) 


are  the  reflection  and  transmission  matrices,  respectively,  for  the 
i'th  layer  referred  to  an  origin  in  the  layer,  and  the  matrices 
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?l  -  exp(±  i  k“  di/2)<S_K  ^  - (2.33) 

g  g"g 

represent  the  propagation  factors  which  account  for  phase  change  and 
attenuation  of  $^(r)  between  the  interlayer  origin. 

Assume  that  scatterings  are  predominantly  in  the  forward  direc¬ 
tion,  and  that  the  absorption  process  is  strong  enough  to  attenuate 
the  wave  amplitude  significantly  within  a  few  layers  from  the  surface, 
then  the  renormalized  forward  scattering  (RFS)  perturbation  theory^ 
can  be  used.  The  amplitude  vectors  A^  can  be  expressed  in  a  series 
of  perturbative  terms  of  order  n, 

A^  -  Z  At(a)  and  A'  -  ZaT^ - (2.34) 

i  n  l  l—i  n  i—i 


A  set  of  recurrence  relation  for  the  wave  amplitudes  between  layers 
can  be  obtained  from  Eq.  (2.31) 


At(n) 

*  P,  T  •  PJ  . 

a+(!° 

+  P+  R4-  P*  A”(n~l) 

i 

i  i  i-1 

i-l 

i  i  i  i 

A-(f 

-  P"  R-4"  P+ 

A+(n> 

+  p”  T  p~  a 

i-1 

i-1  i  i-1 

i-1 

i-1  i  i  i 

,  A-<°) 
Ai 

■0  by  definition. 

and  Ao"^  represents 

- (2.35) 


inside  the  surface  barrier.  If  the  series  expansions  of  these 
equations  are  convergent  at  some  layer  j  for  each  order  (n) , 


is  0  and  then  the  exact  solutions  of  the  transfer  equation 
(2.31)  can  be  obtained.  The  final  backscattered  wave  outside  the 
surface  layer  is  given  by 


(2.36) 


The  inner  most  layer  of  each  pass,  j,  is  chosen  by  the  following 
criterion  of  attenuation  of  the  wavefield 

a  )2cra)  <  0.001  - (2.37) 

I  js 

where  an  is  a  normalization  factor  for  successive  passes 


Convergence  over  successive  passes  is  governed  by  the  same  parameter, 
the  process  stops  at  N  when 


E(a‘In)  )2  )2)<  0.001  - (2.39) 

8  8  g  g 

This  scheme  of  calculation  is  fast  and  efficient,  but  may  break  down 
at  low  energies  where  absorption  is  small  and  back  2.  attering  of  the 
wavefield  is  significant.  The  RFS  method  is  generally  successful  for 
well-separated  layers  with  energies  greater  than  lOev  above  the 


muffin-tin  zero  level. 


Effect  of  thermal  vibration  was  not  considered  in  the  above 
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formulation.  The  thermal  motion  of  the  atom  excites  lattice  waves 
and  distorts  the  regular  arrangement  of  the  atoms,  and  consequently 
reduces  the  intensities  of  the  beams.  The  intensities  of  the  beams 
change  in  the  variation  with  the  energy  and  the  angle  of  the  incident 
beam.  It  is  therefore  an  important  ingredient  in  a  diffraction  cal- 
culation  for  crystals  with  Debye  temperature  considerably  lower  than 
room  temperature.  Including  this  effect  in  a  dynamical  theory  as 
required  in  LEED  is  very  difficult  and  is  presently  done  by  crude 
methods.  Since  the  atoms  of  the  crystal  change  their  positions  by 
an  extremely  small  amount  during  the  time  in  which  an  electron  moves 
into  the  crystal  and  out  again,  it  is  convenient  to  use  the  approxi¬ 
mation  that  the  lattice  motion  is  actually  static  during  the  inter¬ 
action.  We  need  only  consider  the  scattering  of  an  electron  by  a 
crystal  in  which  each  atom  is  displaced  from  its  site  by  a  small 
random  vector  determined  by  the  equilibrium  phonon  distribution. 

We  assume  that  the  electrons  are  scattered  in  the  interiors  of  the 
atoms  rather  than  in  relatively  smooth  potential  towards  the  outside 
and  between  them.  It  is  then  possible  to  calculate  the  scattering 
matrix  t  of  an  atom  on  its  lattice  site  and  use  the  result  to  des¬ 
cribe  the  scattering  of  the  atom  displaced  off  the  site  by  using  the 
same  matrix  t  in  a  shift  sec  of  coordinates.  A  method  introduced  by 
Duke  and  Laramore^,  presently  used  in  most  LEED  calculations,  is  to 
replace  each  vibrating  atom  of  the  crystal  by  a  "blurred  atom".  The 
"blurred  atom"  has  a  scattering  matrix  which  consists  of  the  scat¬ 
tering  matrix  of  the  atom  averaged  over  the  displacements  it  under¬ 
goes  in  the  vibration  motion  of  the  lattice.  The  resulting  lattice 


of  the  blurred  atoms  has  the  full  symmetry  of  the  rigid  lattice  so 
that  the  multiple  scattering  can  be  applied  immediately. 


if  the  scattering  matrix  in  momentum  space  of  a  typical  atom  is 
denoted  by  then  the  scattering  of  this  atom  displaced  by  the 

vector  Ap  is  given  by 


exp ( i (k '  -k)  *Ap) - (2.40) 

and  the  scattering  matrix  of  the  corresponding  blurred  atom  is  given 
by  the  thermal  average  of  the  scattering  matrix,  e.g., 

<t+v£>  av  -  t£v£<exp(i(k"-k)  *Ap)>av - (2.41) 

The  average  occurring  here  is  a  typical  Debye-Waller  factor  and  can 
be  evaluated  to  give 

<t£v*>  av  ■  c£"k  exP<Kk"*“^)  *Ap>  av - (2.42) 

If  the  atom  is  spherically  asymmetric  and  isotropic,  these  relations 

can  be  used  to  define  a  temperature  dependent  "phase  shifts"  5 ^ (T) 

at  temperature  T  for  a  blurred  atom  in  terms  of  the  phase  shift  of 
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5^  of  the  fixed  atom  by 

exp(i6^(T) sin(5^(T) )  i  exp(-2ak2)j  (-2iak2)exp(i5^)x 

sin(5 , )  4tt(2L+1)  (2Z"+1)  /  (21+1)  B2"  (L0,£0) - (2.43) 


with  k  *  2mE/ft  -2Vo 


The  quantity  a  is  given  by 


u 


a  ■  3h 


where  m  is  the  mass  of  the  atom,  kg  is  the  Boltzmann  constant,  and 
is  the  Debye  temperature. 


The  temperature  corrections  can  be  incorporated  into  the  LEED 
calculation  by  substituting  S^by  5^ (T>  in  Eq.  (2.29).  Complete 
theoretical  treatment  of  temperature  effect  in  LEED  intensity  calcu¬ 
lation  can  be  found  in  the  literatures* 


The  computer  block  programs  used  for  our  dynamical  formulation 
of  LEED  intensity  calculation  are  shown  in  Fig.  2.6.  The  details  of 
the  computer  input  and  output  will  be  discussed  in  Chapter  Four. 
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Fig.  2.6  Block  diagram  of  a  dynamic  computer  program 


2.2  Theoretical  Calculation  of  Surface  Scattering  Theory 

2.2.1  Introduction 

When  an  electromagnetic  wave  is  incident  on  a  perfect  plane 
interface  between  two  media,  the  reflected  wave  is  a  function  of  its 
wavelength,  tb«  angle  of  incident  and  the  electrical  properties 
(permittivity,  permeability  and  conductivity)  of  the  two  adjoining 
media . 

The  general  and  exact  solution  of  a  wave  reflected  from  an  ir¬ 
regular  boundary,  the  random  or  periodic  variation  with  respect  to 
height  measured  from  a  certain  mean  plane  is  yet  unknown.  To  find 
the  approximate  solutions  of  such  a  reflected  wave  have  aroused  to 
the  interest  of  physicists  and  engineers  for  the  past  70  years. 

The  foremost  barrier  in  integrated  optics  is  that  the  optical 
propagation  losses  of  the  III-V  compound  semiconductors  cannot  be 
reduced  to  an  acceptable  level.  One  of  the  most  important  sources 
of  loss  of  a  wave  propagating  along  a  thin-film  waveguide  is  the  sur 
face  scattering  effect  from  the  rough  surface  of  the  waveguide  (or 
any  other  structure) .  In  order  to  evaluate  the  performance  of  a 
realistic  waveguide,  it  is  necessary  to  study  its  behavior  if  depar¬ 
tures  from  the  perfect  geometry  occur. 


r  ■ 

ti 


B 
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2.2.2  Theoretical  Calculation  of  Surface  Scattering  of  a 
Symmetric  Slab  Waveguide 

Theoretical  calculation  of  the  surface  scattering  of  a  symmetric 
slab  waveguide  has  been  reported  by  Marcuse^  based  on  mode  conversion 
to  guided  and  radiation  modes.  The  structure  of  a  symmetric  slab 
waveguide  is  simple,  and  could  be  used  to  study  the  radiation  and 
conversion  properties  of  dielectric  waveguide. 

Light  propagation  in  waveguides  occur  in  normal  modes;  the  eigen- 
solutions  of  the  wave  equation  are  obtained  by  solving  the  boundary 
conditions  of  the  waveguide.  There  are  two  types  of  normal  modes: 
"waveguide"  modes,  in  which  the  waves  are  bounded  within  the  guiding 
material,  and  "radiation"  modes,  in  which  the  waves  are  not  bounded 
to  a  dielectric  slab  and  persist  undiminished  through  out  all  space. 
Together  the  "waveguide"  and  "radiation"  modes  form  a  complete  set  of 
eigensolutions  and  have  nondegenerate  eigenvalues,  and  are,  therefore, 
mutually  orthogonal.  The  eigenvalues  of  the  "waveguide"  modes  are 
discrete,  thereas  those  of  the  "radiation"  modes  are  continuous. 

\ 

The  ideal  slab  waveguide  is  sketched  in  Fig.  2.7  (a).  It  resem¬ 
bles  a  longitudinal  cross-section  through  the  claded  optical  fiber. 

We  simplify  the  discussion  further  by  assuming  there  is  no  variation 
of  either  the  waveguide  geometry  or  the  field  distributions  in  y 
direction.  We  express  this  limitation  on  guide  symmetry  by  rhe 
relation 


3/3y  -  0 


(2.45) 
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This  restriction,  Eq.  (2.45),  allows  us  to  decompose  the  field  of  the 
slab  waveguide  into  the  transverse  electric  (TE)  modes  and  transverse 
magnetic  (TM)  modes.  We  begin  with  the  derivation  of  the  properties 
of  the  TE  modes.  By  definition,  the  electric  field  in  the  z  direc¬ 
tion,  E  ,  is  equal  to  zero.  Using  Eq.  (2.45)  and  maxwell's  equations 
we  see  that  only  the  magnetic  field  in  Che  x  and  z  directions,  Hx  and 
H^,  respectively,  and  the  electric  field  in  the  y  direction  E^  are 


non-zero. 

We  can  express  H  and  H 

X  z 

in  terms  of  E  , 

y 

H  - 

X 

-(i/wu)(3E  /3z)  - 

- (2.46) 

and  H  “ 

(i/wu) (3E  /3x)  - 

- (2.47) 

z 

where  w  is  the  frequency  and  ]i  is  known  as  the  magnetic  permeability. 
The  E  -component  is  obtained  as  a  solution  of  the  reduced  wave  equa¬ 


tion: 


(32Ey/3x2)  +  (32Ey/3y2)  +  n2k?Ey  -  0  - (2.48) 


with  n2  ■  e/to - (2.49) 

and  k  ■  w/i7u7  *  2tt/ X 0  - (2.50) 


where  n  is  index  of  refraction,  t  is  the  dielectric  permittivity, 
to  and  Uo  are  the  permittivity  and  magnetic  permeability  in  free 
space,  respectively,  and  Xo  is  Che  wavelength  in  vacuum.  Including 
the  time  and  z-axis  dependent, 


—(2.51) 


39 


exp(i(wt-3z) ) 


and  using  Eqs.  (2.48)  and  (2.51),  we  obtain  the  following  equation, 

(92Ey/3x2)  +  (n2ko-82)  Ey  -  0 - (2.52) 

The  propagation  constant  3  denotes  the  eigenvalue  of  the  normal 
modes . 

According  to  the  orthogonality  relation*2,  any  arbitrary  field 
distribution  can  be  expressed  in  terms  of  the  orthogonal  modes  of  the 
slab  waveguide.  It  is  then  possible  to  write 

Ey  “  ECVEVy  +  E/q(P)  Ey  (P}  dP - (2.53) 

with  C  -  (6  12  wu0P)/“E  E*  dx 
u  u  0  i<x>  y  uy 

and  q (p)  ■  (B/2wu0P)/*EyE  ^  (p)  dx  - (2.54) 

The  first  term  in  Eq.  (2.53)  sums  over  all  TE  waveguide  modes,  while 
the  second  term  is  the  combination  of  summation  and  integration  which 
extends  over  all  radiation  modes.  The  power  P  was  assumed  to  be  the 
same  for  all  the  modes.  The  actual  amount  of  power  that  each  mode 
contributes  to  the  field  expansion  is  determined  by  the  values  of  the 
expansion  coefficients  (C  and  q(p)). 

''J 


For  the  TE  mode,  the  total  power  carried  by  the  field  is  given 
by  the  expression. 


PT  -  P(gJ  Cu|  2  +  Z/?h(p)  |  2  dp  ) - (2.55) 


Our  treatment  of  slab  waveguide  with  imperfect  walls  is  based 
upon  finding  solutions  of  Che  reduced  wave  equation  (2.48)  throughout 
all  space  by  satisfying  certain  boundary  conditions  at  infinity.  It 
is  also  possible  to  express  the  field  of  an  imperfect  waveguide  in 
terms  of  the  modes  of  an  ideal  guide.  We  assume  that  the  walls  of 
the  imperfect  waveguide  obey  the  restriction  of  Eq.  (2.45).  This 
limits  the  boundary  imperfections  to  changes  in  slab  thickness  and 
does  not  allow  variations  in  the  y  direction.  The  departures  of  Che 
waveguide  boundary  have  to  be  very  small  so  that  they  can  be  treated 
by  perturbation  theory.  A  slab  waveguide  is  shown  in  Fig.  2.7  (b) . 
The  distribution  of  the  refractive  index  for  waveguide  with  imperfect 
walls  can  be  indicated  as  follows; 

n2(x,z)«no  (x,z)  +  n(x,z) - (2.56) 

where  n0  describes  the  distribution  of  the  refractive  indices  for 


the  perfect  slab  waveguide. 


g.  2.7  Cross  section  view  of  a  slab  waveguide  (a)  with  perfect 
boundary,  (b)  with  imperfect  boundary. 
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The  values  of  n  are  given  by  the  following  equation. 


for  d+f(z)  <  x 
for  d  <  x  <  d  +  f(z) 

for  -d  +  h(z)  <  x  <d  - (2.58) 

for  -d  <  x  <  -d  +  h(z) 
for  x  <  -d 
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Where  f(z)  and  h(z)  are  wall  distortion  functions.  The  difference 
2  2 

n^-^  may  nor.  necessarily  be  small  for  perturbation  theory  to  be 
applicable.  It  is  sufficient  that  the  region  over  which  n  is  non-zero 
is  narrow. 


C 


We  substitute  Eqs.  (2.53)  and  (2.56)  into  Eq.  (2.48)  and  obtain 


1 


5((32cu/az2)-2isuOcu/3z)4cuk?n)Euy  +  Zf7( 

3  2q(p)/3z2)-2i8(3q(p)/3z)+q(p)k,2n)Ey(p)do-0 - (2.59) 


4 


by  using  the  fact  that  the  modes  E^  and  E^(p)  satisfy  the  wave 
equation  for  the  ideal  slab  waveguide. 


4  We  multiply  Eq.  (2.59)  kv  g  uE^/2w)ja  (u  is  an  integer)  and  inte 

grate  over  the  entire  err ss-section  of  the  guide.  Using  the 
orthogonality  relation,  we  obtain 

4 

02C  /3Z2)  -  2ig  (SC  /3z)  +  I  C  F  (z) 

U  U  U  U  u  'JU 


4 


K 


1 


+Z(°q(p)  Gu(p,z)dp  -  0 - (2.60) 

where  F^z)  =  (8ljlto/2wUop) i^E*uyn  ^x»z)dx - (2-6l> 

and  Gu(p)  -  (3uk5/2wu0P)/Vy  n  (x,z)Ey  (p)dx— (2.62) 
Similarly,  we  obtain 


(32q(p')/3Z2)  -  210'  (3q(p')/3z)  +  §  C^p'.z) 

+Z/°°q(p)  G(p.p-)  dp-  0 - (2.63) 

0 

where  F  (p^z)-(g1co/2wu0P)/“E*(p')r|(x,z)E  dx - (2.64) 

\J  ■AO  y  y 

and  G(p,0-)-(3-kJ/2wUoP)A*(p')n(x,z)Ew(p)dx - (2.65) 

—oo  y  y 
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► 


i 


The  propagation  constant  of  the  radiation  modes  on  p'(a  distance  in 
the  x  direction)  is  indicated  by  Q' .  It  is  apparent  that  the  pertur¬ 
bation,  q,  couples  the  expansion  coefficients  and  q(p)  to  each 
other. 


Eq.  (2.60)  has  the  solution  of 


4 

C  (z)  -  A  +  B  exp(2i8  z)  +  (l/2i8u)  (-/z<i>  (r)dr 

U  UU  U  0USS 


+  exp  (2iS  z)/z  exp(-2iS  Ob  (OdC  - (2.66) 

U  D  U  U 

I 

with  <p  (C)»-EC  F  (z)  -Z/°°q(p)G  (p,z)  do - (2.67) 

u  u  u  uu  0  u 


4 


It  is  advantageous  to  separate  the  coefficient  Cy  into  two  parts. 
The  first  part 


C+  -  A  -  (1/21S  )A>  (5)  dC 

U  U  U  0  u 


■(2.68) 


contributes  to  waves  traveling  in  the  positive  z  direction.  The 
second  part 


C~*(B  +(l/2i8  )/zexp(-2iB  £)4>  (C)d5)exp(2iB  z) - (2.69) 

VI  VI  U  0  VI  VI  u 


contributes  to  waves  traveling  in  the  negative  z  direction. 


The  coefficients  A  and  B  must  be  determined  from  the  initial 
u  u 

conditions.  We  assume  chat  the  lowest  order  TE  mode  is  entering 
longitudinally  from  a  perfect  slab  waveguide  into  a  waveguide  section 
of  length  L,  which  contains  wall  imperfections.  The  imperfect  piece 
of  waveguide  is  again  followed  by  a  perfect  waveguide  of  infinite 
length.  This  mode  allows  us  to  calculate  the  loss  contribution  of  an 
imperfect  waveguide  section  of  length  L.  We  assume  that  at  z»0  only 
lowest  order  TE  mode  travels  in  the  pos-z  direction.  This  implies 


Co  (0)  -  1 


- (2.70) 


and  Cy  (0)  *  0  for  ui*0 


■(2.71) 


We  postulate  that  at  Z*L  no  modes  travel  in  neg-z  direction 


C  (L)  ■  0  for  all  u 


■(2.72) 


The  initial  conditions  of  Eqs.  (2.71)  and  (2.72)  lead  to 


A  -  0 
u  ou 


(2.73) 


<$ou  is  the  Kronecker  delta  symbol,  which  is  zero  for  ui*0  and  unity 
for  u*0.  The  condition  of  Eq.  (2.72)  leads  to 

B  -  (-1/20  )/Lexp(-2iB  (£)  d£ - (2.74) 

u  U  Q  u  u 

The  total  power  loss  can  be  expressed  by  the  following  equation 


AP-PT(L)-PT(0)-Pu?0((l-60u)|C^(L)|2+|c‘(0)!2)+ 

2  00 

{  fo  (|q+(p,L)|2  +  |q~  (P,0)|2)  dp - (2.75) 

The  summation  from  one  to  two  in  Eq.  (2.75)  reminds  us  to  add  the 
losses  contributed  by  the  even  and  odd  radiation  modes. 

Using  this  method,  we  obtain  general  formulas  for  the  power  loss 
in  a  symmetric  slab  waveguide  although  further  approximation  is  needed 
for  practical  application. 

18 

In  the  next  section,  we  introduce  a  simple  theory  for  an 


asymmetric  slab  waveguide. 


2.2.3 


Simple  Theory  of  Surface  Scattering  for  an 
Asymmetric  Slab  Waveguide 
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For  an  asymmetric  slab  waveguide,  the  two  surfaces  of  guiding 
film  scatter  light  differently.  Since  it  is  virtually  impossible  to 
measure  the  scattering  losses  of  the  two  surfaces  separately,  all  the 
surface  properties  of  the  guiding  film  are  lumped  into  a  single  para¬ 
meter.  We  use  this  parameter  to  calculate  the  scattering  losses  of 
different  waveguide  modes. 

A  dielectric  asymmetric  slab  waveguide  is  shown  in  Fig.  2.8. 

Let  us  consider  a  space  of  unit  length  in  the  y  direction  and  assume 
(3/3y-0)  (i.e.,  it  is  infinitely  long  in  the  y  direction).  If  a  light 
wave  is  incident  upon  the  upper  surface  of  the  guiding  region  covering 
a  unit  length  of  the  film  in  the  y  direction,  the  plane  wave  has  a 
width  of  cosd  in  the  direction  parallel  to  the  wavefront.  Considering 
a  TE  wave,  the  incident  power  P  carried  by  the  plane  wave  is  given  by^ 

P  -  (l/2)n1Ey/E77iI7cos0 - (2.76) 

where  Uo  and  £o  are  the  permeability  and  permittivity  of  free  space 

respectively,  0  is  the  incident  angle,  n^  is  the  refractive  index  of 

the  film  and  E^  is  the  field  amplitude.  According  to  the  Rayleigh 
19 

criterion  ,  the  specularly  reflected  beam  from  the  upper  film  sur¬ 
face  has  a  power  P ’ ^ 

P'*(l/2)n^E^/£o/LoCOS0exp-(4TT5^2cost3 /^) - (2.77) 
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where  ^  is  the  wavelength  in  the  guiding  region,  and  5^  is  the 
standard  deviation  of  the  upper  film  surface  roughness  (i.e.,  the 
root-mean-square  deviation  from  its  mean  value) .  The  power  loss  due 
to  the  upper  film  surface  scattering  is  P-P".  Similarly,  the  power 
loss  due  to  the  lower  film  surface  scattering  can  be  obtained  by  re¬ 
placing  <5^2  by  <$^q  equation  2.77.  The  total  power  losses,  k? ,  at 
two  surfaces  of  the  guiding  film  is  therefore 

AP  *  (1/2)0. E^/e  /iiflcos0  (l-exp-K2cos0) 

I  y  o 


*(l/2)n1E^K2  cos39/£o/Uo  - (2.78) 

with  K  -  4ttA1(5J2  +  <5j0)1/2 - (2.79) 


We  have  obtained  a  simple  parameter  K  which  contains  the  surface 
properties  of  the  guiding  dielectric.  It  is  a  dimensionless  parameter 
and  quantitatively  compares  surface  roughness  with  the  optical  wave¬ 
length. 

The  dielectric  slab  waveguide  is  a  useful  model  for  more  compli¬ 
cated  waveguide  structures.  Its  simplicity  allows  us  to  study  the 
properties  of  wave  propagation  in  dielectric  waveguide  without  en- 
cumberance  of  tedious  mathematical  expressions.  However,  in  most 
practical  applications  more  complicated  waveguides  are  used: 
integrated  optics  use  rectangular  strips  of  dielectric  material 
that  are  embedded  in  other  dielectrics. 


2.2.4  The  Surface-Scattering  Losses  of  a  Rectangular 
Dielectric  Waveguide 

A  wave  propagating  along  the  z-direction  in  a  rectangular  di¬ 
electric  waveguide  with  dimensions  a  and  b  (Fig.  2.9)  is  formulated 

and  solved  in  this  section.  An  exact  analytical  treatment  of  this 
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problem  is  impossible,  and  therefore,  Marcatili's  approximate, 
analytical  approach  is  followed.  This  approach  works  only  for  modes 
much  greater  than  cutoff.  For  modes  of  this  type,  its  field  is  con¬ 
fined  to  the  region  of  the  core.  Very  little  field  energy  is  carried 
in  the  surrounding  media.  Let  us  simplify  the  analysis  by  ignoring 
the  field  energy  carried  in  the  shaded  regions  of  Fig.  2.9. 

Maxwell's  equations 

7x  H  *  e0n2  (3E/3t)  and  Vx  E  »  Uo(3H/3t) 

can  be  solved  to  express  the  transverse  field  in  terms  of  longitudinal 
components  Ez  and  Hz>  if  a  traveling  wave  propagating  in  the  Z-direc¬ 
tion  with  incident  angle  9  is  assumed  (Fig.  2.4). 


E  -  (-i/K2)(S(3E  /3x)  +  u0w(3H  /3y)) - (2.80) 

X  Z  2 

H  -  (-i/K2) (S(3H  /3y)  +  we0n2(3E  /3x))  - (2.81) 

y  z  z 

E  -  (-i/K2)(8(3Ez/3y)  -  wUo(3H2/3x))  - (2.82) 


*  ,  r* 


H  -  (-i/K2)(g(3H  /3x)  -  we0n2(3E  /3y)) 
x  z  z 


■(2.83) 


with  K  ■  /(n2k2  -  02) 


where  n  is  refractive  index  of  the  dielectric,  k  is  the  wave  vector 
in  free  space,  w  is  the  angular  frequency  of  the  wave  and  8  is  the 
propagation  constant  in  z-direction.  The  longitudinal  components 
and  Hz  must  satisfy  the  reduced  wave  equation 


(32j>/3x2)  +  (32i|;/3y2)  +  K2^  -  0 - (2.85) 


Two  different  types  of  modes  can  be  supported  in  the  waveguide 


E^  :  polarized  in  x  direction,  and 


E^  :  polarized  in  y  direction 


where  p  and  q  are  positive  integers  Indicating  the  number  of  extrema 

of  the  electric  or  magnetic  field  in  the  x  and  y  direction.  In  the 

E*  modes,  the  main  transverse  field  components  are  E  and  H  .  Within 
pq  x  y 

the  waveguide  each  component  varies  sinusoidally  both  along  the  x  and 
y  directions.  In  region  I,  with  refractive  index  n,  as  shown  in 
Fig.  2.9,  the  wave  equation  can  be  solved  to  give 


E_(x,y)  *  A  cosk  (x+C)  Cosk  (y+n ') - (2.86) 

Z  X  y 


1 


Fig.  2.9  Cross  section  of  a  rectangular  dielectric  waveguide  with 

dimensions  a  and  b.  n^ ,  02 »  n_,  n^  and  n^  are  the  indices 
of  refraction  in  regions  1  to  5,  respectively. 


where  £  and  rTare  the  phase  factors  of  E  with  H  *0,  we  can  obtain 

Z  X 

the  expression  for  E^  from  Eq.  (2.83) 

2  . 

H  "-ACn./B)  /(£o  /Uo)  (k  /k  )  sink  (x+£)  sink  (y+n'’) - (2.87) 

A  y  x  x  y 

Other  components  of  the  wavefunction  can  be  obtained  by  substituting 
Eqs.  (2.86)  and  (2.87)  into  Eqs.  (2.80),  (2.81)  and  (2.82). 


E  »(i  A/k  8) (n2k2-k2)  sink  (x+£)  cosk  (y+n') - (2.88) 

«  xix  x  y 

Ey*  -iA(ky/6)coskx(x+£)  sink^(y+n') - (2.89) 

H  »lA(eo/Uo)n2(k/k  ) sink  (x+£) cosk  (y+n) - (2.90) 

x  i  x  x  y 

with  k2-w2e0Uo  and  n2k2-82«k2  +  k2 - (2.91) 

i-  x  y 


For  well-guided  modes  B^n^k.  This  relation  together  with  Eq.  (2.91) 
result  the  inequalities 

K  <  <  B  and  k  «3  - (2.92) 

x  y 

Owing  to  the  inequalities  in  Eq.  (2.92),  it  is  clear  that  E  «  E  , 

y  z' 

and  Ey  can  be  neglected.  Thus,  there  are  only  four  nonvanishing 


components . 


Similarly,  the  wavefunctions  in  the  surrounding  regions  2,3,4, 
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and  5  with  their  respective  refractive  indices  n2>  n^,  n^,  and  can 

be  obtained.  When  boundary  conditions  are  matched  between  region  1 
and  neighboring  regions  (i.e.,  regions  2,3,4  and  5)  four  transdental 


equations  are  obtained: 

tank  a  *  n2k  (n2s,  +  n2s_) / (n2n2k2  -n*s,s,) - (2.93) 

x  1x35  53  35x  135 

tank^i  »  (n^/n^2  (k^/s^) - -(2.94) 

tank  b  *  k  (s-  +  s,)/(k2  -s_s ) - (2.95) 

y  y  L  4  y  l  4 

tankyn  *  ”s4/,ky  ~  (2.96) 

with  s^  ■  ((n2  -n2  )k2  -k2)  for  i»2,3,4,  and  5  - (2.97) 


The  computer  program  in  appendix,  was  developed  to  solve  k  and  k 

x  y 

by  means  of  Newton's  method.  Once  k  and  k  are  determined,  the 
J  x  y 

propagation  constant  can  be  obtained  from  Eq.  (2.91).  The  total 
power  P  carried  by  mode  in  the  z  direction  given  by 

P  ■  0^3100/(2^0 /£ o )  (//E^dxdy)  -  n^sinQA2/ (4/u0/£q) 

(b+l/s2+l/s^) (a+  (n2n2/s3) (k2+s2) /(n^k2  +  n*s2) 

+  (n2n2/s5)(k2  +  s2)/(n2k2  +njs2))  - (2.98) 


If  a  plane  wave  is  incident  on  the  interface  between  region  1 
and  region  2,  with  an  incident  angle  of  9  covering  the  complete  length 
of  core  in  the  y  direction,  the  plane  wave  has  a  width  of  cos0  in  the 


direction  parallel  to  the  wavefront.  The  incident  power  P  carried 
by  the  plane  wave  is  given  by 

P-  (l/2)n1A2b/e0/ui'  cos0 - (2.99) 

where  A  is  the  field  amplitude.  The  power  loss  AP  at  the  two  inter¬ 
faces  is  given  by 

AP  -(1/2)  n1A2/e0/Uocos9(l-exp(-K2cos29))b 
-(l/2)n^A2K2  cos30  b^So/Ua - (2.100) 

with  K  -  4tt/(622  -h5J3)1/2 

where  6i2  and  are  c^e  standard  deviation  of  the  interfaces  between 
media  1  and  2,  and  between  media  1  and  3,  respectively.  Dividing 
Eq.  (2.100)  by  Eq.  (2.98),  yields  the  power  attenuation  per  unit  length 
of  the  waveguide 

a-  K2(cos39/sin9)  (a+  (n2n2/s.j)((k2+s2) /(n^k2  +n^s2^ 

+  (n2n2/ Sj) (k2+s2) / (n^k2  +njs2))"1+(b+i/s2+l/s4)*1 - (2.101) 

This  is  the  loss  caused  by  surface  scattering.  This  loss  is  expressed 
as  a  product  of  four  independent  factors.  The  first  factor  K,  depends 
solely  on  the  surface  properties  of  interfaces  and  is  dimensionless. 

The  second  factor  involving  9  depends  on  the  waveguide  mode.  The 
third  and  fourth  factors  are  the  effective  waveguide  thickness  and 


width,  respectively. 


CHAPTER  3 
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EXPERIMENTAL  METHODS 


3.1  Auger  Electron  Spectroscopy  (AES) 

Auger  electron  spectroscopy  is  a  nondestructive  method  of  sur¬ 
face  chemical  analysis.  It  explores  the  electronic  energy  levels  in 
atoms  and  solids  by  analyzing  the  energies  of  the  electrons  ejected 
from  a  specimen  bombarded  with  an  electron  beam.  Some  of  the  ejected 
electrons  result  from  the  "Auger  process".  The  "Auger  process" 
involves  any  electron  de-excitation  in  which  this  de-excitation 
energy  is  transferred  to  and  causes  the  ejection  of  a  second  electron, 
the  "Auger  electron".  Because  of  the  discrete  nature  of  most  elec¬ 
tronic  energy  levels,  the  energy  of  the  Auger  electron  i3  character¬ 
istic  of  the  atom  from  which  it  was  ejected.  If  the  Auger  electron 
is  created  within  a  few  atomic  layers  of  the  specimen  surface,  it 
can  migrate  to  the  surface  and  escape  into  the  vacuum  without  loss 
of  energy.  Thus,  the  Auger  electrons  produce  peaks  in  the  secondary 
electron  energy  distribution  and  the  measurement  of  these  energies 
Indicate  the  type  of  atoms  present. 

21 

Auger  electrons  were  first  discovered  in  1925  by  Pierre  Auger 

These  electrons  are  produced  by  the  "Auger  process"  as  shown  in 
22 

Fig.  3.1  for  a  silicon  atom.  The  electronic  energy  levels  of  the 
single  ionized  atoms  are  listed  on  the  left,  taking  the  Fermi  level 
Ep,  as  zero  energy.  The  x-ray  levels,  K,  L,  M. . . ,  are  also  labeled 


on  che  rlghc.  Fig.  3.1  shows  the  atom  immediately  after  a  K-shell 
ionization  by  a  primary  electron.  The  Auger  process  initiates  when 
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an  outer  electron  fills  this  "hole",  as  shown  for  an  electron,  and 
the  energy  released  by  this  transition  is  either  emitted  as  a  photon 
or  given  to  another  electron.  If  this  energy  is  sufficient,  an  elec¬ 
tron  can  be  ejected  from  the  atom,  as  illustrated  for  another  electron 
at  1*2  2*  The  process  just  described  is  called  the  KL^L2  3  Auger 
transition.  The  energy  E  of  the  ejected  electron  can  be  approximated 
by  the  diagram  of  Fig.  3.1.  The  energy  released  is  E^-E^;  however, 
the  ejected  electron  must  expend  the  energy  (E£2  3+  <P)  to  escape  the 
atom.  <j>  is  the  work  function  of  the  material  and  E^2  3  /E^2  3^s  c*ue 
to  the  extra  positive  charge  of  the  atom.  Because  of  this  positive 
charge,  E£2  3  should  be  approximately  equal  to  the  I*2  3  ionization 
energy  of  the  next  heavier  element,  or  E£2  ^  3  >  wber®  z 

is  the  atomic  number  and  A- 1  accounts  for  the  extra  charge.  We  then 


have 
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E(Z)-Ek(Z)-EL1(2)-EL2  3(Z+A)-$ - (3.1) 

Eq.  (3.1)  can  be  generalized  to  any  Auger  event  involving  the  levels 
WXY.  The  actual  measured  energy  will  have  an  additional  term 
(-<p.H>)  which  is  the  difference  between  the  work  functions  of  the 

A 

energy  analyzer  $  arid  the  material  under  investigation.  This  gives 

A 

Ewxy(Z)-Ew(Z)-Ex(Z)-Ey(Z+A)-^A - (3.2) 


In  practical  applications,  Auger  peaks  can  be  greatly  enhanced 
by  recording  the  derivative  of  the  secondary  electron  energy  distri¬ 
bution.  This  technique  improves  the  sensitivity  and  resolution  of 
the  measurement  and  removes  the  large  slowly  changing  background  upon 
which  the  Auger  signal  is  superimposed.  Thus,  the  conventional  Auger 
spectrum  is  the  function  dN(E)/dE.  The  locations  of  the  Auger  peaks 
are  specified  by  the  energies  of  the  negative  excursions  of  the  peaks 
in  the  Auger  spectrum.  The  peak-to-peak  magnitude  of  an  Auger  peak 
in  a  differential  spectrum  generally  is  directly  related  to  the  sur¬ 
face  concentration  of  the  element  which  produces  the  Auger  electrons. 
Quantitative  analysis  may  be  accomplished  by  comparing  the  peak 
heights  obtained  from  an  unknown  specimen  with  those  from  pure  ele¬ 
mental  standards  or  from  compounds  of  known  composition. 

Since  the  escape  depth  without  energy  loss  for  Auger  electrons 
in  the  practical  energy  range  is  typically  less  than  a  few  atomic 
layers.  Auger  electron  spectroscopy  (AES)  is  primarily  a  surface 
analytical  technique.  The  sensitivity  of  the  Auger  wechnique  is 
determined  by  the  transition  probability  of  the  Auger  transitions 
involved,  and  the  incident  beam  current  and  energy.  The  limit  of 
detection  is  around  1 Z  of  a  monolayer,  and  the  sensitivity  varies 
less  than  a  factor  of  10  over  all  elements  except  hydrogen  and 
helium  which  cannot  be  detected. 

Our  experimental  arrangement  for  obtaining  the  Auger  spectra  is 
shown  in  Fig.  3.2.  Several  specimens  were  mounted  on  the  manipulator 
and  sequentially  rotated  into  the  analysis  position.  The  analysis 


was  carried  out  with  the  excitation  beam  normal  to  the  specimen  sur¬ 
face.  The  basic  function  of  each  of  the  electronic  units  is  indicated 
in  Fig.  3.2.  The  Electron  Gun  Control  operates  the  electron  gun,  the 
Auger  System  Control  supplies  the  modulated  scan  voltage  for  the 
cylindrical  analyzer,  and  the  lock-in  Amplifier  provides  phase  sensi¬ 
tive  detection  of  the  electron  multiplier  output.  The  gain  of  the 
electron  multiplier  can  be  adjusted  for  optimum  performance  using  the 
Electron  Multiplier  Supply.  Readout  of  the  Auger  spectrum  is  provided 
by  an  oscilloscope  or  X-Y  recorder.  Auger  data  was  taken  in  the 
dN(E)/dE  mode.  Electronic  differentiation  was  accomplished  with  a 
velocity  analyzer  by  superimposing  a  small  a.c.  voltage  on  the  energy 
selecting  voltage  and  synchronously  detecting  the  output  of  the  elec¬ 
tron  multiplier.  The  operating  conditions  were  2V  peak-to-peak  a.c. 
modulation,  an  incident  electron  current  of  10-lSpA,  and  an  incident 

electron  energy  of  5KeV,  unless  otherwise  specified.  The  system  is 

.9 

capable  of  continuously  maintaining  pressures  less  than  1x10  torr. 

Ion-sputtering  was  accomplished  by  using  an  ion  gun  as  shown  in 
Fig.  3.2  which  produced  an  Ar  ion  beam  in  a  static  Ar  pressure  of 
5x10  ^  torr.  The  ion  and  electron  beams  were  mutually  perpendicular 
and  were  approximately  aligned  to  focus  on  the  same  area  of  the 
testing  surface.  Ion-sputtering  can  be  used  along  for  cleaning  pur¬ 
poses  or,  when  combined  with  AES,  for  analyzing  the  surface  to  obtain 
the  chemical  depth  profiles. 


3.2  Laser  Scattering 
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The  laser  scattering  technique  ,  shown  in  Fig.  3.3,  was  per¬ 
formed  to  characterize  the  surface  geometry  (i.e.,  degree  of  "smooth¬ 
ness")  of  GaAs  substrates.  An  average  "Bidirectional  Reflectance 

Function"  (BRDF)  number  was  used  to  describe  the  degree  of  surface 

24 

smoothness  of  the  tested  sample.  The  BRDF  number  is  defined  as 

BRDF»(Pri/Pi)  (cos0isin4)i)  (AcosQ^sin^/R^) - (3.3) 

where  P^  is  the  scattered  power  collected  by  a  photometer  with  aper¬ 
ture  area  A,  is  the  power  Incident  on  the  sample,  and  is  the 
incident  beam  angle  measured  from  vertical  axis  Z  (as  shown  in  Fig. 
3.3).  <j>r  denotes  the  scattered  beam  angle  measured  from  Z  axis,  0^ 
denotes  the  incident  beam  angle  measured  from  horizontal  axis  Y  and 

©r  is  the  angle  of  the  reflected  beam.  The  diameter  of  the  incident 

2 

laser  beam  was  0.5mm  ,  as  shown  in  Fig.  3.4.  The  sample  holder  is 

automatically  controlled  to  have  translation  and  rotation  movements 

in  the  X-Z  plane  which  allows  the  incident  beam  to  scan  over  an  area 
2 

of  9mm  around  the  center  of  each  sample.  Five  locations  of  each 
sample,  depicted  in  Fig.  3.4,  were  first  rotationally  scanned  by  the 
incident  beam.  The  BRDF  values  vs.  rotational  angle  were  plotted, 
and  from  these  curves  the  directions  at  which  the  maximum  and  minimum 
scatter  occurred  was  determined.  Each  sample  was  then  scanned  trans- 
lationally  along  Che  maximum  and  minimum  directions  separately,  and 
then  the  average  BRDF  value  was  obtained.  The  average  BRDF  value 
quantitatively  characterizes  the  degree  of  surface  roughness. 


The  parameter  K  discussed  in  section  2.2.4  can  be  characterized  by 
the  average  BRDF  value.  , 


Fig.  3.3  A  schematic  diagram  of  laser  scattering  measurement. 
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3 . 3  Anodic  Oxidation 

25 

An  Anodic  oxidation  scheme  is  suggested  for  use  in  forming  a 

rib  waveguide.  The  waveguide,  which  consists  of  a  lightly-doped 
15  -3 

(10  cm  )  GaAs  substrate  with  a  rib  etched  into  the  lightly-doped 
GaAs  layer,  can  be  used  to  investigate  the  surface-scattering  losses 
of  III-V  semiconductor  waveguides. 

The  experimental  setup  employed  for  anodic  oxidation  is  shown 
in  Fig.  3.5.  A  quartz  or  pyrex  beaker  with  a  platinum  cathode  was 
used  to  contain  an  electrolyte  of  D. I. -water  with  PH  adjusted  to  2.5 
by  the  addition  of  phosphoric  acid.  The  GaAs  samples  are  either 
partly  or  fully  immersed  in  the  electrolyte  while  held  by  an  aluminum 
clamp.  Either  a  constant  current  or  a  constant  voltage  supply  can  be 
selected  as  the  power  source. 

The  growth  of  a  surface  oxide  on  GaAs  is  carried  out  by  polar¬ 
izing  the  GaAs  sample  positively  with  respect  to  the  platinum  elec¬ 
trode.  Both  the  oxide  thickness  and  GaAs  consumed  in  the  oxide 

o 

growth  are  proportional  to  the  anodizatloc  voltage  (20  and  12  A/V, 

26 

respectively)  .  The  sample  can  then  be  etched  in  KOH  as  the  surface 
oxide  is  cracked  away  to  give  the  desired  rib  shape. 


Fig.  3.5  Experimental  setup  for  anodic  oxidation. 
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3.4  Electrical  Measurements 

3.4.1  Resistivity  Measurement 

27  28 

Resistivity  measurements  carried  out  by  the  four-point  probe  ’ 

29 

and  Van-der  Pauw  methods  were  used  to  determine  the  doping  concen¬ 
trations  of  the  GaAs  substrates  as  stated  in  the  previous  section. 

The  results  obtained  by  both  methods  showed  good  agreement. 

The  four-point  probe  method,  as  shown  in  Fig.  3.6,  is  the  most 
common  method  for  measuring  resistivity.  A  small  current  (I)  (1  to 
100  mA  in  our  case)  supplied  from  a  constant  current  source  is  passed 
through  the  outer  two  probes  and  the  voltage  (V)  between  the  inner 
two  probes  is  measured.  The  spacing  S  between  the  probe  is  one  milli¬ 
meter.  For  a  thin  sample  with  thickness  W  (in  the  present  case  the 

-2  10 
thickness  is  5x10  cm) ,  the  resistivity  is  given  by 

P-(V/I)  *CF*W  ohm-cm - (3.4) 

where  CF  is  the  correction  factor  shown  in  Fig.  3.5.  The  doping 
concentration  (N)  for  an  n-type  substrate  is  obtained  by 

N«l/(u  qp)  cm - (3.5) 

U 

where  y  is  the  electron  mobility,  and  q  is  the  electron  charge. 


The  Van-der  Pauw  method  was  carried  out  by  making  four  small 
contacts  along  the  circumf erence  of  the  substrate  as  shown  in  Fig. 
3.7  and  by  measuring  the  two  resistance  R^  qj  and  R^  ^  and  the 
thickness  W  of  the  sample.  The  resistance  R,_  __  is  defined  as  the 

Ad  ,  LU 

voltage  V  across  the  contacts  D  and  C  per  unit  current  through  the 
contacts  A  and  B.  The  resistivity  is  given  by 

0-2.72^^^  +  R^^)  f(RAB,CD/SBC,DA)  (3*6) 


where  the  function  f(R^  CD^^BC  DA^  is  s^own  itl  3.7. 


CORRECTION  FACTOR  CF 


Fig.  3.6  Correction  factor  for  resistivity  measurement 
using  a  four-point  probe. 


RAB,CD/R8C,DA 
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Fig.  3.7  Resistivity  measurements  using  Van-Der  Pauw  method,  (a)  a 
sample  of  arbitrary  shape  with  four  contacts  at  vary  place 
(b)  the  function  f  used  for  determining  the  resistivity. 


3.4.2 


Current-Voltage  (I-V)  Characteristics 


The  circuit  diagram,  shown  in  Fig.  3.8,  was  used  to  measure  the 
I-V  characteristic  under  forward  bias  for  GaP  Schottky  diodes.  The 
sample  was  mounted  in  a  shielded  box  with  a  spring  loaded  copper  wire 
making  contact  to  the  top  electrode.  The  current  through  the  diode 
was  measured  using  an  electrometer  (PAR  135) .  The  current  ranged 
from  fractions  of  a  nanoampere  to  microamperes  for  GaP  with  aluminum 
or  gold  Schottky  contacts.  The  I-V  curves  were  then  plotted,  and  the 
slope  and  Y-axis  intercept  of  this  curve  was  used  to  determine  the 
barrier  height  (<$,  )  and  the  diode  factor  (n) . 


3.4.3  Capacitance-Voltage  (C-V)  Measurement 

The  C-V  characteristics  under  reverse  bias  for  GaP  Schottky 
diodes  were  measured  using  a  Boonton  capacitance  meter  operating  at 
1  MHz  as  shown  in  Fig.  3.9.  The  d.c.  bias  was  applied  using  a  regu¬ 
lated  power  supply  connected  to  the  back  of  the  instrument.  The  a.c. 
signal  superimposed  on  d.c.  bias  was  applied  to  the  GaP  Schottky 
diode. 

30 

For  a  Schottky  barrier  ,  at  reverse  bias  and  high  frequency, 
the  capacitance  per  unit  area  C  is  given  by 

C-(qesNQ/2(V-<0bn))1/2 - (3.7) 

where  q  is  the  electron  charge,  is  dielectric  constant  of  GaP, 

is  the  carrier  concentration  in  the  GaP  substrate  and  is  the 

2 

barrier  height.  The  plot  of  1/C  vs.  voltage  V  should  be  a  straight 
line  with  slope  of  l/qN-C,.,  and  having  an  intercept  on  the  X-axis 
equal  to  4>bQ. 


C— V  measurement  using  Boonton  Capacitance  Meter 


3.4.4  Capacitance-Frequency  (C-f)  and 

Conductance-Frequency  (G-f)  Measurements 

C-f  and  G-f  measurements  for  GaP  Schottky  diodes  at  zero  bias 
were  performed  using  a  lock-in  amplifier  (ITHACO  model  339)  along 
with  a  photometric  preamplifier  (Par  model  184)  as  shown  in  Fig.  3.10. 
In  the  experimental  setup,  a  capacitor  C  was  used  to  prevent  the  bias 
voltage  supply  from  interfering  with  the  reference  signal  of  the 
lock-in  amplifier.  The  resistance  R  provides  a  high  impedance  path 
to  the  lock-in  amplifier's  reference  signal  and  acts  as  a  current 
limiting  resistor  in  case  the  sample  short  circuited,  thus  preventing 
overloading  the  lock-in  amplifier.  The  RC  network  acting  as  a  high 
pass  filter  does  not  reduce  the  useable  range  of  frequencies  (1.5Hz- 
150KHz) .  The  lock-in  amplifier  was  calibrated  using  a  suitable 
standard  capacitor  read  at  full  scale  deflection.  The  magnitude  of 
the  reference  signal  was  20  millivolts.  The  sample  current  is  re¬ 
solved  into  in-phase  and  quadrature  components,  and  a  d.c.  voltage 
proportional  to  the  current  provided  at  the  output  of  the  lock-in 
amplifier.  The  voltage  corresponding  to  the  in-phase  and  quadrature 
components  represented  the  conductance  G  and  the  Capacitance  C, 
respectively. 


C-f  and  G-f  measurements  using  lock-in  amplifier. 
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CHAPTER  4 

ATOMIC  STRUCTURE  OF  GaP  (110)  FACE 

4 . 1  Introduction 

/ 

III-V  compound  semiconductors  are  used  widely  in  solid  state 
electronic  devices  and  in  chemical  processes  which  include  adsorption 
and  catalysis.  Typical  examples  of  III-V  semiconductor  devices  which 
cannot  be  realized  using  elemental  semiconductors  include  light  emit¬ 
ting  diodes  (LED's),  Gunn  effect  oscillators,  and  high  frequency 
field  effect  transistors.  Recognizing  these  applications,  it  becomes 
necessary  to  obtain  a  quantitative  understanding  of  the  free  surface 
properties  of  compound  materials  so  that  we  may  develop  a  fabrication 
technology  of  the  same  quality  as  that  which  exists  for  silicon-based 
devices . 

It  is  common  knowledge  that  low  energy  electron  diffraction 

(LEED)  is  a  tool  of  great  potential  in  the  study  of  surface  structures 

and  can  be  used  quantitatively  to  determine  the  atomic  structure  of 

single  crystal  surfaces.  The  surface  atomic  structure  of  GaAs 

(HO)^-^  and  InP  (110)  has  been  successfully  characterized  by 

LEED.  The  atomic  geometry  of  the  GaAs  (110)  face  was  used  as  the 

starting  point  to  study  its  surface  chemisorption,  oxidation  and 
36—38 

passivation  ,  since  an  accurate  knowledge  of  surface  atomic 


geometry  is  essential  to  most  studies  of  surface  related  properties. 
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GaP,  with  an  energy  gap  of  2.3  eV,  is  one  of  the  III-V  compounds 

39 

used  for  light  emitting  diodes  (LED's)  ,  and  for  highly  efficient 

40-43 

negative  electron  affinity  (NEA)  electron  emitters  .  It  has  been 

found  that  the  atomic  structures  and  surface  conditions  of  GaP  are 

41  42 

important  to  the  adsorption  of  Cs  on  the  GaP  (100)  ,  GaP  (110)  , 

43 

and  GaP  (111)  crystal  faces.  In  this  chapter,  the  results  of  an 

atomic  surface  structure  determination  for  the  (110)  face  of  GaP  is 
44 

reported  .  A  quantitative  determination  of  the  structure  for  the 
GaP  (110)  face  is  obtained  by  comparing  a  dynamic  computer  aided 
simulation  of  LEED  intensities  with  those  measured  experimentally. 

The  LEED  patterns  and  experimental  results  for  the  GaP  (111)  face 
are  analyzed  and  reported,  and  the  surface  structure  of  the  GaP  (110) 
face  is  compared  with  that  of  the  GaAs  (110)  face. 


4.2  Experimental  Considerations 

All  experimental  measurements  of  the  LEED  intensities  were  per¬ 
formed  in  a  standard  ultrahigh  vacuum  LEED/AES  system.  The  schematic 
diagram  of  the  system  is  shown  in  Fig.  4.1.  It  consists  of  a  four- 
grid  hemispherical  LEED  optics  with  fluorescent  screen  display,  a  5-KV 
cylindrical  mirror  analyzer  and  a  rotatable,  temperature-controlled, 
manipulator.  Three  sets  of  Helmholtz  coils  were  used  to  align  the 
incident  electron  beam.  LEED  intensities,  recorded  from  the  screen 
with  a  telescopic  spot  photometer  and  plotted  out  on  an  X-Y  recorder, 
were  subsequently  corrected  for  variation  in  primary  beam  current  and 
background  intensity.  The  background  intensity  of  the  screen  was 
recorded  and  subtracted  from  the  beam  intensities.  The  beam  inten¬ 
sities  were  then  normalized  to  the  LEED-gun  electron  emission  current. 
The  diffraction  patterns  were  recorded  photographically.  Intensity- 
voltage  profiles  were  taken  at  normal  incidence,  which  was  assured 
by  determining  that  symmetrically  equivalent  beams  provided  the  same 
intensity  profiles  to  within  a  few  percent.  Since  this  procedure  is 
only  accurate  to  within  a  few  degrees,  sample  alignment  is  a  major 
source  of  uncertainty  in  the  intensity  data. 

The  GaP  sample  was  mounted  on  the  sample  holder  of  the  manipu¬ 
lator.  After  insertion  into  the  ultrahigh  vacuum  (UHV)  chamber,  the 
surface  wa3  cleaned  by  1  Kev  sputtering  using  argon  ions  for  20  min 
and  was  thermally  annealed  at  550° C  for  1  hour.  A  sharp  lxl  LEED 
pattern  was  obtained  for  the  GaP  (110)  face  as  shown  in  Fig. 4. 2(d). 


79 

An  Auger  profile  showed  that  both  oxygen  and  carbon  (the  principle 
contaminants)  peaks  were  reduced  to  the  limit  of  the  analyzer  resolu¬ 
tion  indicating  that  the  surface  was  clean  and  well  ordered. 

LEED  intensities  were  recorded  at  both  room  temperature  and  low 
temperatures.  The  sample  was  kept  at  90K  during  data  recording  by 
flowing  liquid  nitrogen  through  the  sample  holder,  and  both  room 
temperature  and  low  temperature  data  were  used  to  analyze  the  surface 
structure. 

Only  one  set  of  GaP  (110)  LEED  experimental  intensity-voltage 
45 

(I-V)  curves  were  available  at  the  time  of  our  data  acquisition. 

In  our  work  a  much  larger  experimental  data  base  was  obtained,  and 
was  in  general  agreement  with  that  recorded  previously.  LEED  profiles 
for  a  total  of  16  inequivalent  nonspecular  beams  were  recorded  and 
analyzed.  The  reproducibility  of  the  I-V  curves  was  satisfactory. 

Fig.  4.2  (a),  (b)  and  (c)  show  the  notation  for  the  LEED  beams.  The 
arrangement  of  the  surface  atoms  of  GaP  (110)  face  viewed  from  the 
3ide  and  top  positions  are  shown  in  4.2  (a)  and  (b) ,  respectively, 
and  the  notation  used  in  labelling  the  (h,k)  beam  is  shown  in  (c) . 

The  LEED  pattern  taken  at  164eV  with  liquid  nitrogen  cooling  the 
sample  (95K)  is  shown  (d) . 


Turbo 

Pump 


Fig.  4.1 


Schematic  diagram  of  the  ultra  high  vacuum  (UHV) 
system  equipped  with  LEED,  AES  and  ion-sputtering 
gun  for  surface  analysis. 
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Fig.  4.2 
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(a)  Side  view  of  the  GaP  (110)  face.  (b) 

Top  view  of  surface  atoms  for  the  undistorted 
GaP  (110)  face.  (•)  Ga,  (o)  P,  d  is  the 
nearest  neighbor  distance,  (c)  Tne  notation  used 
in  labeling  the  (h,k)  beam,  (d)  LEED  pattern  of 
GaP  (110)  face  taken  at  low  temperature. 


4.3  Model  and  Theoretical  Considerations 

An  approximate  multiple-scattering  model  of  the  diffraction 

process,  described  in  section  2.5,  was  used  to  construct  the  computer 

program  for  calculating  the  theoretical  1-7  curves.  In  the  intensity 

program,  the  intralayer  multiple  scattering  was  treated  exactly  using 

46 

the  direct  summation  method  ,  while  for  the  interlayer  multiple 
scattering,  an  approximate  fast  scheme  called  renormalized-forward- 
scattering  (RFS)^’^  perturbation  theory,  was  used.  In  the  atomic 
potential  program,  the  electron-ion  core  interaction  is  described  by 
a  one-electron  muffin-tin  potential.  To  calculate  the  muffin-tin 

potential  for  atoms  in  a  crystal  lattice,  according  to  the  description 

49  50  51 

of  Matheiss  ’  ,  we  started  with  the  atomic  wavefunction3  and  then 

included  a  correction  for  ionic  materials  (i.e.,  Ga+P  )  by  a  madelung 

52 

type  summation  in  real  and  reciprocal  space.  Slater  exchange  was 

used  for  the  contribution  of  exchange  interaction.  The  muffin-tin 
+  -  0 

radii  of  Ga  and  P  were  chosen  to  be  1.182A.  This  is  the  half  of 
the  nearest-neighbc  :  distance  of  GaP.  The  phase  shifts  are  then  ob¬ 
tained  by  radial  integration^  of  the  Schrodinger  equation  for  the 
separated  atomic  potential.  The  first  five  phase  shifts  (£*0,1-4) 
were  used  in  the  dynamical  LEED  intensity  calculations  through 
Eq.  (2.1-36).  The  electron-electron  interaction  was  incorporated 
via  a  complex  "optical  potential".  The  inner  potential  (Vor)  was 
taken  to  be  8  eV.  The  inelastic  damping  factor  (imaginary  potential) 


was  set  equal  to  5  eV. 


The  consequence  of  the  thermal  lattice  vibrations  (electron- 
phonon  interaction)  was  included  in  the  LEED  intensity  calculations 
by  correlating  the  input  phase  shifts  with  the  temperature  dependent 
phase  shifts  through  Eq.  (2.1-57).  A  Debye  temperature  of  319K  was 
used  for  GaP.  Experimental  LEED  intensities  were  very  Insensitive  to 
the  temperature  factor.  Both  room  temperature  and  low  temperature 
(95K)  LEED  intensity  curves  of  the  GaP  (110)  face  were  recorded  and 
essentially  identical.  This  indicates  that  the  approximate  Debye 
model  will  not  attribute  any  error  to  the  determination  of  the  GaP 
(110)  surface  structure. 


4.4  Surface  Structure  of  GaP  (110)  Face 

The  arrangement  of  the  surface  atoms  on  the  ideal  GaP  (110) 

surface  in  the  direct  lattice  is  the  same  as  that  of  GaAs  except  that 

the  As  atoms  are  replaced  by  P  atoms.  The  dimensions  of  the  surface 
0  °54 

unit  mesh,  a«3.864A  and  b-5.4505A  as  shown  in  Fig.  4.2  (b) ,  are  very 
close  to  those  of  GaAs.  Since  the  GaP  crystal  is  quantitatively  simi¬ 
lar  to  GaAs  with  respect  to  its  crystal  structure  and  chemical  pro¬ 
perties,  it  is  reasonable  to  assume  that  the  GaP  (110)  face  has  the 
same  surface  reconstruction  as  that  of  the  GaAs  (110)  face.  No  frac¬ 
tional  order  spots  were  observed  as  shown  in  Fig.  4.2  (d)  and  the 
ratio  of  the  dimensions  of  the  surface  unit  mesh  is  1.41.  Agreement 
with  the  ideal  bulk  value  of  1.414  is  within  1 7.  which  suggests  that 
the  surface  unit  mesh  is  undistorted  with  respect  to  that  in  the 
bulk  unit  cell.  The  twofold  symmetry  along  the  b  axis  (Fig. 4. 1(b)) 
in  the  direct  lattice  of  the  (110)  face  is  directly  observed  through 
the  symmetry  of  the  intensity  profiles  of  the  (±h,k)  beams.  This 
suggests  that  a  relative  horizontal  displacement  of  the  group  III  and 
the  group  V  atoms  in  the  surface  layer  is  unlikely.  The  intensities 
of  the  (1,0)  and  (1,0)  beams  are  expected  to  be  extremely  weak  when 

they  are  diffracted  from  a  bulk  structure  because  the  (±1,0)  beams 

32 

are  "quasi-forbidden  beams"  .  This  can  be  explained  by  the  formula 
for  the  scattering  amplitude 

Fhk"  ijafj  ^»9^exP(27ri^hxij  +k^ij  +U+c°s0hk)Zj A)) - (4.1) 


microcopy  resolution  test  chart 
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where  is  the  atomic  form  factor  of  the  atom  in  the  j'th  layer, 

and  y^  are  the  coordinates  of  the  i'th  atom  in  the  j'th  layer, 

is  the  distance  from  the  surface  to  the  j'th  layer,  8^  denotes  the 

scattered  angle  of  (h,k)  beam  and  is  the  Incident  electron  wave  length. 

For  the  (1,0)  and  (1,0)  beams  of  GaP  (110),  h«±l,  k*0  and  z»id.  The 
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summation  in  Eq.  (4.1)  then  gives: 

F±i,o“(fGa‘fP)/(1‘exp(1+cos9hk)2lTidA) - (4*2) 

Equation  (4.2)  shows  that  scattering  from  rows  of  Ga  and  rows  of  P 
atoms  is  out  of  phase  and  should  cancel.  Since  the  form  factor  of 

Qa  p 

Ga  and  P  (f  and  f  )  should  be  similar,  Eq.  (4.2)  predicts  that  the 
intensities  of  the  (1,0)  and  (1,0)  beams  should  be  weak  if  there  is 
no  displacement  in  the  position  of  the  atoms  on  the  (110)  face.  The 
observed  intensities  of  (±1,0)  beams  are  very  strong  as  shown  in  the 
experimental  LEED  patterns  (Fig.  4.2(d)),  and  in  the  intensity  pro¬ 
files  (Fig.  4.3).  The  discrepancy  between  theoretical  and  experi¬ 
mental  results  suggests  that  there  is  a  rearrangement  in  the  positions 
of  the  atoms  on  the  (110)  surface.  The  observed  unreconstructed  LEED 
patterns  and  the  excellent  symmetry  of  the  (±h,k)  beams  rules  out  the 
possibility  of  surface  unit  mesh  reconstruction  and  the  horizontal 
displacement  of  Ga  and  P  atoms  with  respect  to  each  ocher.  A  verti¬ 
cal  contraction  of  the  top  surface  layer  or  a  vertical  displacement 
of  Ga  and  P  atoms  with  respect  to  each  other,  which  will  leave  the 
surface  unit  mesh  vnchanged,  is  the  only  remaining  possibility  for 
surface  rearrangement.  All  of  the  experimental  results  obtained 
from  the  GaP  (110)  face  are  similar  to  those  obtained  from  the  GaAs 
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(110)  face 
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Based  on  Che  information  of  GaAs  (110)  face"3  ’  ’  ,  we  assumed 
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the  bond-rotation  model  with  w*271'  (Fig.  4.1(a))  was  the  best  candi¬ 
date  for  surface  atomic  geometry  of  the  GaP  (110)  face.  This  assump¬ 
tion  turned  out  to  be  incorrect,  as  will  be  shown  later.  In  the  study, 
other  possible  surface  structures  analyzed  in  detail  include  (a)  pure 
bond-rotation  with  w»20°,27°,  and  34. 8°,  (b)  pure  top  layer  spacing 
contraction,  Ad,  and  (c)  a  bond-rotation,  w,  with  a  top  layer  spacing 
contraction,  Ad.  A  total  number  of  17  conceivable  surface  structures 
were  analyzed  in  detail.  These  structures  are  summarized  in  Table  4.1. 

I-V  curves  obtained  for  the  first  16  inequivalent  nonspecular 
beams  emerging  from  the  sample  were  compared  with  the  theoretical 
curves  calculated  from  the  17  different  structures  as  listed  in 
Table  4.1.  Most  of  the  17  proposed  structures  yielded  very  poor 
results  when  we  compared  the  theoretical  curves  with  the  experimental 
curves.  Therefore,  these  proposed  structures  could  not  be  considered 
as  the  real  surface  structure.  Only  3  out  of  17  structures  provided 
relatively  good  agreement  with  the  experimental  curves,  and  these 
structure  models  were  (a)  a  bond-rotation  angle,  w»27° ,  with  a  5Z 

9 

top  surface  layer  spacing  contraction,  d*-0.097A,  (b)  a  pure  bond- 

o 

rotation  angle,  w»27° ,  plus  Ad»0,  and  (c)  a  pure  0.097A  contraction 
of  top  layer  spacing  with  w«0.  These  theoretical  curves  and  experi¬ 
mental  curves  are  compared  as  shown  in  Figs.  4.3  to  4.6  for  the  (0,1), 
(1,0),  (0,1)  and  (1,1)  beams  respectively  for  the  first-order  beams. 

The  comparison  between  the  recorded  and  calculated  curves  for  the 


high-order  beams  such  as  (0,2),  (1,2)  and  (2,1)  beams  are  similar  to 
the  first-order  beams.  Obviously,  the  curves  (c)  in  Figs.  4.3  to  4.6 
are  not  in  good  agreement  with  the  experimental  curves.  Thus,  the 
model  of  a  pure  5%  top  layer  spacing  contraction  could  not  be  consi¬ 
dered  for  real  surface  structure.  Both  curves  (a)  and  (b)  are  only 
in  fair  agreement  with  the  experimental  curves.  From  the  peak  posi¬ 
tions  and  the  line  shape  of  the  (1,0)  beam  in  Fig.  4.4,  the  (0,1) 
beam  in  Fig.  4.5,  and  the  (1,1)  beam  in  Fig.  4.6,  model  (a)  is  favored 
rather  than  model  (b)  as  the  real  surface  atomic  structure  of  GaP 
(110)  face. 

The  generally  satisfactory  agreement  of  most  of  the  16  beams  in 
magnitude,  peak  positions,  and  line  shape  between  the  recorded  and 
calculated  curves  of  model  (a)  led  to  the  conclusion  that  P  atoms  are 
rotated  outward  from  the  GaP  (110)  plane  and  Ga  atoms  are  rotated  in¬ 
ward  from  this  plane  by  an  angle  of  27°.  Furthermore,  the  surface 

layer  is  compressed  by  5%  such  that  the  top  layer  spacing  is  reduced 

o 

by  approximately  0.1A. 

No  quantitative  analysis  or  the  GaP  (111)  surface  atomic  struc¬ 
ture  has  been  reported.  A  lxl  LEED  pattern  was  observed  by  Miyao, 
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et  al,  and  then  a  2x2  reconstruction  pattern  was  reported  by  Van 
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Bommel  and  Crombeen  .  In  our  study,  a  clear  lxl  threefold  LEED 
pattern  of  GaP  (111)  was  first  observed  in  a  UHV  chamber  with  back- 

_Q 

ground  pressure  of  1x10  torr  as  shown  in  Fig.  4.7,  and  a  2x2 
reconstruction  pattern  was  then  obtained  in  a  better  background 
pressure  of  4xl0~*°  torr.  The  extra  (2x2)  reconstructed  beams  were 
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weak  and  disappeared  within  hours  after  the  sputtering  and  annealing. 
Further  study  of  this  phenomenon  is  needed.  The  dynamical  calculation 
and  analysis  is  in  progress,  but  complete  results  will  not  be  reported 
here. 


One  interesting  experimental  result  is  worth  mentioning.  A  LEED 

pattern  indicating  diffraction  from  facets  was  observed  on  the  GaP 

(111)  face  after  annealing  the  crystal  to  600°C.  The  facet  structure 

has  been  reported  previously  and  it  was  concluded  that  the  facet 

pattern  is  due  to  reflections  from  three  equivalent  (110)  planes^. 

The  facet  spots  faded  away  after  heating  to  650° C  and  spots  of  the 

fundamental  lxl  pattern  brightened.  The  same  results  were  reported 
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by  Miyao,  et  al.  .  No  explanation  of  these  results  was  offered  by 
them.  We  attempt  to  explain  this  phenomenon  as  follows.  A  large 
amount  of  P  atoms  start  to  evaporate  from  the  (111)  face  and  leave 
an  almost  pure  Ga  layer  on  the  surface  when  the  sample  is  heated  to 
650°C.  At  this  temperature,  the  Ga  layer  is  a  melted  metallic  layer 
and  crystallizes  along  the  GaP  (111)  structure  during  the  crystal 
cooling  period.  An  epitaxial  Ga  layer  is  thus  grown  on  GaP  (111) 
surface.  The  ratio  of  the  Ga/P  Auger  peak  increase  from  0.9  to  5.05 
after  heating  to  650°C  which  gives  the  support  to  our  explanation. 
This  interesting  result  indicates  that  it  should  be  relatively  easy 
to  obtain  a  good  Schottky  contact  on  the  GaP  (111)  face. 


Table  4.1 


Surface  structure  for  GaP  (110)  face. 


Ad  is  the  percentage  of  the  top  layer  spacing 
contraction  and  w  is  the  bond  rotation  angle 
with  the  Ga  species  rotated  inward  and  P 


with  Ad-0 


GaP(1lO) 

(0,1) 

*4m-OS7l 

•• 

/ 

0 

— y 

(b) 

Ad*  0 

a 

u>-  27° 

/ 

■J  (a! 

Ad— j097A- 

INTENSITY  (arbitrary  units) 


GaP  (110) 

(1,0) 

(e) 

Ad»-.097A 

s 

W3 

(b) 

Ad*0 

94 

4.5  Discussion  and  Conclusion 

The  LEED  pattern  of  the  GaP  (110)  face  starts  to  fade  away  after 

-9 

very  few  hours  in  a  1x10  torr  UHV  chamber  while  the  GaAs  (110)  LEED 
pattern  can  stay  sharp  for  days  in  the  same  vacuum  environment.  In 
both  cases,  no  reconstruction  LEED  spots  were  observed  during  the 
period  that  they  remained  in  the  chamber.  This  implies  that  the 
chemisorption  of  the  residual  gas  molecules  on  both  the  GaP  (110)  and 
the  GaAs  (110)  faces  is  a  random  process.  An  oxygen  Auger  peak  was 
detected  on  the  former  after  the  sample  remained  in  the  chamber  for 
hours,  while  it  took  days  for  the  latter  to  pick  up  any  detectable 
amount  of  oxygen  in  the  same  vacuum  environment.  These  indicate  that 
the  GaP  (110)  face  is  more  reactive  and  less  stable  than  the  GaAs 
(110)  face. 

No  attempt  has  been  made  to  investigate  the  efficiency  of  the 

Zanazzl  and  Jona  r  factors'^  as  a  quantitative  guide  to  the  trend  in 

the  theoretical  intensities.  It  is  very  useful  to  use  r  factors  in 

refining  and  quantifying  an  already  good  theory-experiment  comparison. 

It  is  less  useful  to  use  these  r  factors  for  the  first  qualitative 

model  selection.  More  extensive  data  and  analysis  including  r  factors 

could  lead  to  further  refinements  of  the  structure  model.  In  a  re- 
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cent  report  by  C.  B.  Duke,  et  al.  ,  they  utilized  r  factors  to  assess 
the  analysis.  The  resulting  best-fit  structures  consist  of  single 

layer  reconstruction  characterized  by  a  bond-rotation  angle  of  w*259 

o 

±3°  and  a  0.1±0.05  A  contraction  of  the  topmost  layer  spacing  which 
is  similar  to  our  proposed  model.  Our  result  turns  out  to  be  a  good 


first  approximation. 
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In  conclusion,  our  LEED  study  has  indicated  that  the  atomic 

structure  of  the  GaP  (110)  face  is  quite  similar  to  that  of  GaAs 

(110)  regarding  the  rippled  geometry.  On  both  the  GaP  (110)  and  GaAs 

(110)  surfaces,  the  group  V  atoms  are  rotated  outward  from  the  surface 

plane  and  Ga  atoms  are  rotated  inward  from  this  plane  by  an  angle  of 

27°  in  order  to  preserve  the  bond  lengths  between  the  nearest-neighbor. 

a 

The  topmost  layer  spacing  of  the  GaP  (110)  face  is  reduced  by  0.1A. 

The  stability  study  of  LEED  patterns  concludes  that  the  GaP  (110) 
face  is  less  stable  than  the  GaAs  (110)  face.  The  complete  surface 
structure  analysis  of  the  GaP  (111)  face  is  premature  and  needs  further 
study.  It  is  suggested  that  a  good  Schottky  contact  should  be  rela¬ 
tively  easy  to  grow  on  the  GaP  (111)  face. 


CHAPTER  5 


STRUCTURE  MODEL  OF  III-V  COMPOUND  SEMICONDUCTOR  SCHOTTKY  BARRIERS 

5.1  Introduction 

Many  investigations  have  been  done  on  metal-semiconductor  inter¬ 
faces,  and  numerous  theoretical  approaches  have  been  used  to  explain 

the  basic  physical  mechanism  for  the  metal-semiconductor  interaction 
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and  formation  of  the  Schottky  barriers  .  Modem  surface  analyti¬ 
cal  techniques  have  been  applied  to  characterize  the  metal-semicon- 
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ductor  interfaces  and  the  results  strongly  suggest  that  the 

physical  structure  of  a  real  Schottky  barrier  is  more  complicated 
than  that  of  a  perfect  planar  junction  between  the  metal  and  semi¬ 
conductor.  Defects  at  or  near  the  surface  of  the  semiconductor  are 
the  important  mechanism  in  determining  the  barrier  height  and  metal- 
semiconductor  interaction.  In  many  cases,  this  can  dominate  the 
electrical  properties  of  Schottky  barrier.  Interface  compound  for¬ 
mation  and  the  metal-induced  surface  states  are  observed^  at  the 
representative  metal-semiconductor  interfaces  and  it  is  concluded 
that  local  charge  redistribution  rather  than  any  intrinsic  surface 
states  of  the  semiconductor  determine  Schottky  barrier  formation  at 
these  interfaces.  The  "defect"  model  ’  ^  and  the  "unified"  model  * 
of  all  Schottky  barriers  have  been  proposed  to  explain  the  Fermi-level 
pinning  and  other  phenomena  of  a  real  Schottky  barrier.  Although  the 
Schottky  diode  has  important  applications  in  practical  devices,  and 
the  formation  of  Schottky  barriers  has  been  studied  for  the  past 


several  decades,  both  che  physical  and  electrical  properties  of  this 
device  are  not  completely  understood. 
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In  this  chapter,  a  macroscopic  structural  model  of  the  real 
Schottky  barrier  is  suggested.  Based  on  the  structural  model,  an 
electrical  equivalent  model  is  derived  to  describe  the  fundamental 
behaviors  of  a  real  Schottky  diode,  including  the  origin  of  the  Fermi 
level  pinning,  the  barrier  heights,  and  the  deviation  of  the  diode 
factor  from  its  ideal  value. 

As  suggested  at  the  conclusion  of  Chapter  4,  the  GaP  (111)  face 
was  chosen  to  be  the  semiconductor  surface  for  Schottky  barrier 
studies.  Schottky  barriers  were  fabricated  by  evaporating  metal  on 
a  GaP  (111)  face  to  demonstrate  the  structural  model.  Capacitance- 
voltage  (C-V) ,  capacitance-frequency  (C-f) ,  current-voltage  (I-V) 
and  conductance-frequency  (G-f)  measurements  were  then  performed  to 
analyze  the  electrical  properties  of  the  diodes.  Auger  electron 
spectroscopy  (AES)  studies  were  used  to  correlate  the  chemical  com¬ 
position  of  the  metal-semiconductor  interface  to  the  structure. 
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5.2  Structural  Model 

Metal-induced  surface  states  and  interface  compound  formation 

at  metal-semiconductor  interfaces  were  observed*^.  These  phenomena 

were  associated  with  microscopic  dipoles  at  the  intimate  contacts 

which  accounted  for  the  macroscopic  Schottky  barrier  heights.  Modern 

surface  analytical  techniques  have  been  applied  to  study  the  metal- 
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semiconductor  interfaces  and  the  results  ’  showed  a  strong  dif¬ 
fusion  and  inter-atomic  diffusion  between  the  metal  and  semiconductor. 
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It  was  suggested  that  the  driving  force  for  dissolution  of  the  III-V 

compound  and  its  movement  into  the  metal  is  the  heat  of  condensation 

of  the  metal  on  the  surface  of  semiconductor.  The  energy  transferred 

from  the  metal  atoms  during  the  condensation  must  be  dissipated  by 

the  semiconductor.  If  this  energy  is  comparable  to  or  considerably 

larger  chan  the  heat  of  formation  of  the  semiconductor,  the  momentarily 

excited  atoms  at  the  semiconductor  surface  will  gain  enough  energy  to 

break  their  bonds  and  move  into  the  metal.  This  atomic  interdiffusion 

process  could  create  a  thin  amorphous  layer  which  is  comprised  of  the 

metal  and  semiconductor  atoms  between  the  metal  and  semiconductor. 

A  detailed  discussion  of  the  condensation  mechanism  can  be  found  in 
69 

the  literature 


Strong  interdiffusion  of  the  metal  and  semiconductor  atoms  was 
observed  in  the  AES  studies  of  Au-GaP  and  Al-GaP  Schottky  barriers. 
Based  on  the  experimental  results,  a  metal-amorphous  film-semicon¬ 
ductor  (MAS)  structure  as  shown  in  Fig.  5.1  (a)  is  suggested  to  be 
the  physical  structural  model  of  a  real  Schottky  barrier.  The  thick- 


ness  of  the  amorphous  film,  W,  is  in  the  range  of  50A  to  500A  depending 
upon  fabrication  details,  and  the  amorphous  film  is  a  random  mixture 
of  metal  and  semiconductor  atoms.  The  a.c.  equivalent  circuit  of  the 
physical  structural  model  under  reverse  bias  is  shown  in  Fig.  5.1  (b) . 

C  and  R  are  the  capacitance  and  the  resistance  of  the  amorphous 

A  A 

layer,  respectively,  is  the  depletion  capacitance  of  the  semi¬ 
conductor,  and  Rg  is  the  series  resistance  of  Schottky  diode.  The 
d.c.  equivalent  circuit  of  the  model  under  forward  bias  is  shown  in 
Fig.  5.1  (c) .  At  low  current  injection,  the  I-V  characteristics  of 
an  amorphous  semiconductor  is  a  hyperbolic  sine  function.  At  high 
current  injection,  the  hyperbolic  function  can  be  approximated  as; 

I«IoAexp(VA/nAVT) - (5.1) 

where  VT  -  kT/q 

where  Iq.  is  a  constant,  V  is  the  applied  voltage  across  the  amor- 

A  A 

phous  film,  k  is  Boltzmann's  constant,  T  is  the  temperature  in  degrees 
Kelvin,  q  is  the  magnitude  of  electron  charge,  and  nA  is  a  parameter 
characterizing  the  properties  of  the  film,  including  its  thickness. 
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Fig.  5.1  (a)  is  the  structural  model  of  a  real  Schottky  diode. 

(b)  is  the  a.c.  equivalent  circuit  of  the  real 
Schottky  diode  under  reverse  bias. 

(c)  is  the  d.c.  equivalent  circuit  of  (a)  under 
forward  bias. 


5.3  OuCline  of  Experiments 

Schottky  diodes  described  In  this  chapter  were  prepared  on  single 
crystal  n-type  GaP  (111)  faces  doped  with  Te(N^  -  3x10^).  The  GaP 
crystals  were  1  in. -diameter  by  0.03  in. -thick  wafers,  and  polished 
on  both  (111)  and  (111)  faces.  Eight  samples  were  cut  from  the  same 
wafer.  Low  energy  electron  diffraction  (LEED)  studies  were  conducted 
on  the  samples  before  device  fabrication  to  ensure  that  ohmic  contacts 
were  deposited  on  (ITT)  face.  Ohmic  contacts  were  attached  first  by 
alloying  Au  Si  ^  Co  the  (111)  face  of  the  crystal.  The  alloying 
was  done  at  600°C  in  a  nitrogen  gas  atmosphere  for  three  minutes. 

Low  contact  resistances,  (<1Q) ,  were  obtained  in  most  of  the  diodes. 

Prior  to  the  evaporation  of  rectifying  contacts,  the  surface  of 
each  sample  was  etched  in  a  solution  of  H^SO^tH^O^iD.I.-water  (3:1:1). 
The  ohmic  contacts  were  protected  with  Apiezon  wax  during  this  pro¬ 
cess.  After  removing  the  Apiezon  wax  with  trichloroethylene  and 
rinsing  in  acetone,  methanol  and  then  D. I. -water,  the  substrates  were 
blown  dry  with  nitrogen,  and  then  placed  in  the  vac-ion  vacuum  system 
immediately. 

Rectifying  contacts  were  deposited  on  the  GaP  (111)  face  by  a 

vacuum  evaporation  of  Al  or  Au  from  resistance  sources  with  a  back- 
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ground  pressure  in  the  low  10  torr  range.  The  temperature  of  the 
samples  was  retained  at  250° C  throughout  the  evaporation.  This  ele¬ 
vated  substrate  temperature  reduced  the  possibility  of  interlayer 
formation  between  the  metal  and  semiconductor  due  to  surface  contam- 
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ination.  The  evaporation  sources  were  covered  by  a  shutter  during 
the  cleaning  preheating  period,  and  the  shutter  was  opened  to  allow 
for  a  metal  deposition.  A  deposition  period  of  10  seconds  was  used 
for  all  the  samples  except  for  sample  No.  3,  which  had  a  considerably 
longer  evaporation  time.  A  much  thicker  amorphous  interlayer  was 
obtained  in  sample  No.  3  due  to  the  prolonged  evaporation  cycle,  as 
will  be  noted  later  in  the  AES  profile. 

The  geometries  of  the  rectifying  contacts  were  controlled  by 
0.05  inch  thick  molybdenum  masks,  and  circular  contacts,  0.03  inch  in 
diameter,  were  used  for  C-V,  I-V,  and  AES  measurements. 

A  series  of  C-V,  C-f,  G-f  and  I-V  measurements  as  described  in 
section  3.4  were  carried  out  for  at  least  three  diodes  on  each  sample, 
after  the  sample  had  been  removed  from  the  evaporation  station.  All 
the  measurements  were  performed  in  a  shielded  box.  For  the  reverse 
biased  C-V,  C-f  and  G-f  measurements,  a  lock-in  amplifier  was  used. 

C-f  and  G-f  measurements  were  carried  out  at  zero  bias  for  different 
frequencies.  For  the  forward  bias  I-V  measurement,  an  electrometer 
and  a  digital  multimeter  were  used.  The  chemical  depth  profiles  were 
obtained  by  using  AES  and  ion  milling  in  an  UHV  chamber  after  the 
electrical  measurements. 


5.4  Results  and  Discussion 

The  I-V  characteristics  curves  for  samples  No.  3  through  No.  6 
and  sample  No.  8  were  plotted  in  Fig.  5.2.  The  diode  factors  n  were 
derived  from  the  extrapolation  of  the  linear  region  of  these  curves 
and  are  listed  in  Table  5.1.  Significant  deviation  of  the  diode  fac¬ 
tor  n  from  its  ideal  value  of  unity  is  observed  on  all  samples,  which 
indicates  that  the  current  transport  between  Schottky  barriers  is  not 
purely  thermionic.  This  deviation  can  be  attributed  to  the  existence 
of  the  amorphous  film  in  the  diode,  which  behaves  as  a  nonlinear  re¬ 
sistor  at  high  current  injection  as  modeled  in  Fig.  5.1  (c) .  A  more 
detailed  characterization  of  this  film  at  forward  bias  will  be  dis¬ 
cussed  later.  The  series  resistance  of  the  diodes  were  obtained  from 
the  I-V  curves  at  large  forward  bias.  The  diode  capacitance  per  unit 

area,  C,  was  obtained  from  measurements  of  the  capacitance,  conduc- 

2 

tance  and  diode  series  resistance.  1/C  vs.  V  curves  remained  linear 

for  a  large  region  of  bias  voltages  as  shown  in  Fig.  5.3  for  five 

different  samples.  The  values  of  carrier  concentration  were  calcu- 

2 

lated  from  the  slope  of  1/C  vs.  V  in  the  reverse  bias  region.  The 
Schottky  barrier  heights,  of  all  samples  were  obtained  from  both 

C-V  and  I-V  measurements.  These  values  of  are  listed  in  Table 
5.1. 


From  Table  5.1,  we  noted  that  the  barrier  heights  obtained 

by  I-V  measurements  are  generally  consistent  and  agree  with  data  ob- 
74 

tained  by  others  .  The  barrier  heights  obtained  by  C-V  measurements, 


however,  are  always  much  larger  than  those  deduced  from  I-V  methods. 


l/C  -V  curves  of  GaP  (111)  Schottky  diodes  with  various  metal  films. 
Barrier  heights  and  doping  concentration  of  the  diodes  obtained  from 
the  curves  are  listed  in  Table  S.l. 


Table  5.1 


Results  of  C-V  and  I-V  measurements. 


and  this  is  consistent  with  the  results  generally  reported^ 

In  some  cases,  the  barrier  heights  measured  by  the  C-V  method  were 
larger  than  the  energy  bandgap  of  GaP,  (see  Table  5.1),  where  the 
barrier  heights  of  sample  Nos.  3  and  6  are  3.1eV  and  2.35eV,  respec¬ 
tively.  This  behavior  was  also  observed  and  reported  by  Amith,  et 
al^.  It  was  proposed  that  a  thin  native  oxide  or  a  contaminate 
interfacial  layer  existed  between  the  metal  and  the  semiconductor  to 
explain  this  behavior.  It  did  not,  however,  explain  the  behavior 
observed^  in  diodes  fabricated  on  newly  cleaved  surfaces  in  an  UHV 
chamber,  where  the  contamination  or  oxidation  of  the  semiconductor 
was  eliminated.  Brillson,^ ’ ^ ,  reported  observing  chemical  reactions 
and  charge  redistribution  at  metal-semiconductor  junctions,  and  the 
formation  of  reacted  interlayers  with  frequency-dependent  dielectric 
constants  significantly  different  from  their  semiconductor  counter¬ 
parts,  on  Al-CdS  and  Al-CdSe  diodes.  Our  studies  indicated  that  the 
reacted  Interlayer  revealed  many  characteristics  of  an  amorphous 
semiconductor,  and  therefore,  a  metal-amorphous  film-semiconductor 
structural  model  is  suggested.  Our  experimental  results  obtained 
from  real  Schottky  diodes  tend  to  support  this  model  completely,  as 
will  be  described. 

The  AES  depth-composition  profiles  of  sample  Nos.  3,  5,  6  and  8 
are  shown  in  Fig.  5.4  through  5.7,  respectively.  In  all  cases,  no 
detectable  amount  of  contaminants  were  shown  at  the  metal-semicon¬ 
ductor  interfaces.  This  result,  however,  cannot  rule  out  the 
possibility  of  the  existence  of  contamination  one  percent  or  less 
than  one  percent  of  a  monolayer,  which  is  below  the  limit  of  Auger 


SPUTTERING  TIME  (min) 

AES  depth  -  composition  profile  of  sample  No. 3. 
Interdiffusion  region  was  defined  as  the  amorphous  film  (W) 


SPUTTERING  TIME(min) 

AES  depth  -  composition  profile  of  sample  No.  5. 
Interdiffusion  region  (W)  was  defined  as  the  amorphous  film 


SPUTTERING  TIME( 

omposition  profile  of  sample  No.  6 
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signal  detectability.  Fermi-level  pinning  caused  by  such  contamination 
is  possible.  An  interlock  system,  which  has  the  capability  of  fabri¬ 
cating  Schottky  diodes  in  an  UHV  chamber  on  a  newly  cleaved  or  ion- 
sputtering  cleaned  semiconductor  surfaces  to  avoid  such  contamination, 
is  currently  underdeveloped.  A  small  amount  of  oxygen  and  carbon 
contamination  was  detected  at  the  surface  region  of  the  metal  layer, 
but  this  should  play  no  role  in  the  diode  performance.  Strong  atomic 
interdiffusion  was  observed  for  all  samples.  The  Interdiffusion 
region,  W,  the  width  of  the  amorphous  film  between  the  metal  and 
semiconductor  layers,  is  defined  to  start  at  the  point  where  the 
concentrations  of  metal  and  semiconductor  are  both  50%,  and  to  end 
at  the  point  adjacent  to  the  semiconductor  where  the  concentration  of 
semiconductor  is  95%,  as  indicated  in  the  AES  profiles.  It  will  be 
noted  that  W  is  not  well  defined,  but  since  the  amorphous  film  is  a 
mixture  of  metal  and  semiconductor  atoms,  if  the  concentration  of 
metal  atoms  is  greater  than  50%,  metallic  properties  would  be  ex¬ 
pected  to  dominate  the  characteristics  of  the  mixture.  The  end  point 
of  the  amorphous  film  is  defined  on  the  basis  that  a  few  atoms  of 

metal  condensed  on  the  semiconductor  could  randomize  the  surface 
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structure  of  the  semiconductor  and  generate  enough  defect  states 
to  pin  the  Fermi-level. 

o 

The  ion  sputtering  rate,  8A/min,  was  carefully  calibrated  on  a 
known  thickness.  The  sputtering  rate  and  the  interdiffusion  width  W 
were  combined  to  experimentally  estimate  the  thickness  of  the  amor¬ 
phous  film  in  all  the  samples,  and  the  resulting  values  of  W  are 


listed  in  Table  5.2. 


The  barrier  heights,  (<£.  )  obtained  using  the  I-V  method  and  the 

bn 

carrier  concentrations  (N^)  obtained  by  the  C-V  method,  were  used  to 
calculate  the  depletion  capacitance  (C^)  of  the  semiconductor  at  zero 
bias.  The  total  capacitance  Co  and  the  total  resistance  Ro  of  the 
diode  were  obtained  from  C-f  and  G-f  measurements,  respectively.  The 
values  were  used  together  with  Co  and  Ro  to  calculate  the  capaci¬ 
tance  C  and  the  resistance  R  of  the  amorphous  film  utilizing  the 

A  A 

equivalent  circuit  of  the  diode  under  reverse  bias  shown  in  Fig.  5.1 

(b) .  Values  of  C^  and  R^  for  all  of  the  samples  at  1  KHz  are  listed 

in  Table  5.2.  The  dielectric  constants  £  of  the  film  for  all  the 

r 

samples,  listed  in  Table  5.2,  were  estimated  by  values  of  the  capac¬ 
itance  C^  and  the  width  W.  The  average  value  of  the  dielectric 
constant  is  7.23±0.31.  It  is  clear  that  the  dielectric  constant  of 
sample  5  as  listed  in  Table  5.2  is  much  smaller  than  that  of  the 
other  samples.  This  is  probably  due  to  the  error  introduced  in 
estimating  the  film  thickness  of  sample  5,  and  because  of  this  un¬ 
certainty,  the  value  of  dielectric  constant  of  sample  5  was  not 
included  in  the  calculation  of  the  average  dielectric  constant  for 
the  amorphous  film. 

We  found  that  both  R^  and  er  are  frequency-dependent.  Figs.  5.8 
and  5.9  show  the  R  (f)  and  C  (f)  curves  of  various  samples.  The  re- 

A  A 

lation  between  the  amorphous  resistor  R^  and  frequency  is  given  by 

f”°.  The  values  of  m  were  calculated  and  indicated  in  Fig.  5.8. 

The  frequency-dependent  nature  of  RA  and  reveal  the  amorphous 

characteristics^  of  the  interlayer.  The  exponent  m  has  been  de- 
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rived  ’  theoretically  as 


Table  5.2 


Characteristics  of  the  amorphous  film 
found  in  the  GaP  (111)  Schottky  diodes. 


■  Low-frequency  (1  KHz)  capacitance  of  amorphous  film  at  zero  bias 
R  *  Low-frequency  (1  KHz)  resistance  of  amorphous  film  at  zero  bias 

A 

W  ■  Thickness  of  amorphous  film 

*  Dielectric  constant  of  amorphous  film 


Sample 

Number 

Metal 

CA(pf) 

Vn> 

W(A) 

e 

r 

3 

Al(Ni) 

657 

8.38xl05 

504 

7.5 

5 

Al(Au) 

1113 

5.54xl05 

133 

3.3 

8 

Al(Cu) 

1489 

4.32xl05 

217 

7.3 

6 

Au 

1183 

1.02xl06 

259 

6.9 

f 


FREQUENCY  f(Hz) 

amorphous  resistance  verse  frequency  at  zero  bias 


FREQUENCY  UHz) 

Frequency  dependence  of  the  dielectric  constant  for  the  film. 


n“d(lna)  /d(w)»l-4/ln(Vp^/w) - (5.2) 

where  a  is  the  conductivity  of  the  amorphous  semiconductor,  w*2irf  is 
the  frequency  of  the  measuring  signal,  and  Vp^  is  a  phonon  frequency 
(typically  on  the  order  of  10  S  ).  Theoretically,  m  is  calculated 
to  be  0.8  at  approximately  lOKHz.  The  experimental  value  of  m  was 
measured  to  be  0.79±0.22  at  a  frequency  of  lOKHz.  This  result  supports 
the  hypothesis  regarding  the  amorphous  nature  of  the  interlayer. 

The  dielectric  constant  of  amorphous  GaP  was  calculated  by  using 
the  amorphous  film  width,  and  the  capacitance  CA  obtained  by  C-f 
measurement.  This  information  was  not  available  in  the  literature. 

The  frequency  dependent  curves  of  the  dielectric  constant  are  plotted 
in  Fig,  5.9  for  various  samples. 

The  experimental  values  of  diode  factors  reported****’^  for  real 
Schottky  diodes  are  generally  much  larger  than  the  ideal  value  of 
unity.  Our  study  reinforced  these  results  and  the  average  value  of 
diode  factors  was  found  to  be  1.41+0.10.  We  attribute  the  higher 
value  of  the  diode  factor  to  the  existence  of  the  amorphous  film  in 
Che  diode.  At  forward-bias,  the  d.c.  equivalent  circuit  of  a  real 
Schottky  diode  is  shown  in  Fig.  5.1  (c) .  The  total  voltage  drop 
across  the  real  diode,  V,  consists  of  two  components,  VQ,  and  V^. 

is  the  voltage  drop  on  the  amorphous  region  and  is  the  voltage 
drop  on  the  ideal  diode  D.  We  neglect  the  voltage  drop  on  the  con¬ 
tact  resistor  because  it  is  very  small  when  compared  with  and  V  . 


The  total  current  I  in  the  diode  D  can  be  given  by 


I-I0De]",<VVT) - <5'3) 

where  1^  is  the  reverse  saturation  current  of  an  ideal  diode  and  n^ 
is  the  ideal  diode  factor  which  is  unity.  Referring  back  to  section 
5.2,  the  I-V  characteristics  of  the  amorphous  film  can  be  approximated 
by  an  exponential  form  at  high  current  injection  as  given  by  Eq.  (5.1) 
Combining  Eqs.  (5.1)  and  (5.3)  and  with  the  fact  that  V*VD+VA  yields; 

I»I0exp  (V/  VT)  *IQexP  (V/nVT) - (5 . 4) 

"here  I0’I0D(WI0D)(°*/”> 

«• 

where  H-n^-Hi^-l-H^is  the  diode  factor  generally  measured  experimen¬ 
tally.  Eq.  (5.4)  demonstrates  that  a  higher  value  of  diode  factor 
usually  measured  can  be  attributed  to  the  existence  of  an  amorphous 
film. 


The  I-V  curves  of  sample  Nos.  3  and  S,  as  shown  in  Fig.  5.10 

indicate  that  sample  No.  3  needed  an  additional  voltage  (AV)  in  order 

to  obtain  the  same  current  level  as  that  of  sample  S.  As  listed  in 
Table  5.1,  it  is  clear  that  the  barrier  heights  of  sample  Nos.  3  and 
8  obtained  from  I-V  measurements  are  almost  identical.  It  is  reason¬ 
able  to  assume  that  the  voltage  drop  on  the  ideal  diode  V  is  the  same 

for  both  samples.  The  only  difference  between  them  is  that  the  amor- 

o 

phous  film  in  sample  No.  8,  (W*217  A)  is  much  thinner  than  that  of 


sample  No.  3,  (W*504  A).  Sample  No.  3  must  have  a  smaller  amorphous 
resistance  than  that  of  sample  No.  3,  and  the  additional  voltage  drop, 
(AV)  on  the  latter  is  actually  dropped  across  the  amorphous  resistor 
R^.  Therefore,  a  plot  of  AV  against  the  current  density  J  obtained 
from  Fig.  5.10,  shown  in  Fig.  5.11  is  actually  the  indirect  plot  of 
I-V  characteristics  of  the  excessive  amorphous  resistor  in  sample 
No.  3,  as  compared  with  that  of  sample  No.  8.  An  exponential  func¬ 
tion  is  shown  in  Fig.  5.11,  and  this  experimentally  verifies  the  I-V 
characteristics  of  the  interlayer  modeled  by  Eq.  (5.1).  The  above 
argument  is  proof  that  the  film  which  exists  between  the  metal  and 
semiconductor  of  a  real  Schottky  diode  has  characteristics  of  an 
amorphous  nature. 

Based  on  the  MAS  structural  model,  extrinsic  surface  states  on 
the  semiconductor  surface  would  be  created  before  an  amorphous  film 
is  formed.  Our  model  strongly  supports  the  defect  model  which  ex¬ 
plains  the  origin  of  Fermi-level  pinning  of  real  Schottky  barriers, 
and  can  be  combined  with  the  defect  model  to  explain  the  transition 
in  interface  behavior*’*’  between  covalent  and  ionic  semiconductors. 

Fig.  5.12  illustrates  the  significant  of  chemical  reactivity  in 

Schottky  barrier  formation.  S,  indicated  in  Fig.  5.12,  is  the  index 
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of  interface  behavior  as  defined  by  Kurtin,  et  al  .  AH,  in  Fig.  5.12, 

is  the  chemical  heat  of  semiconductor  formation  which  can  emphasize 

the  role  of  interface  reactivity.  From  Fig.  5.12,  it  can  be  noted 

that  covalent  semiconductors  have  a  heat  of  formation  smaller  than 

40Kcal/mole,  whereas  the  heat  of  formation  for  a  metal  on  a  semicon- 
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ductor  is  generally  larger  than  60Kcal/mole  .  When  metallic  atoms 


are  deposited  on  the  semiconductor,  the  energy  (>60Kcal/mole)  trans¬ 
ferred  from  the  atoms  to  the  semiconductor  during  the  condensation  is 
larger  than  the  heat  of  semiconductor  formation  (<40Kcal/mole) .  This 
energy  must  be  dissipated,  and  thus,  the  bonds  of  the  semiconductor 
surface  atoms  would  break.  The  result  would  be  a  random  interdiffusion 
of  metal  and  semiconductor  at  the  interface  and  the  formation  of  amor¬ 
phous  interlayer,  as  well  as  creating  surface  defects  on  the  semi¬ 
conductor.  The  Fermi-level  is  then  pinned  by  these  surface  defects. 

On  the  other  hand,  semiconductors  with  a  heat  of  formation  larger 
than  60Kcal/mole  are  much  more  stable,  because  the  condensation  heat 
of  a  metal  is  generally  not  large  enough  to  break  the  atomic  bonds 
of  these  semiconductors  during  the  Schottky  barrier  formation.  The 
formation  of  an  amorphous  layer  and  creation  of  surface  defects  are 
not  likely  to  happen  in  these  semiconductors.  The  magnitude  of  the 
barrier  heights  formed  on  these  semiconductors  is  highly  dependent 
upon  the  difference  between  the  work  function  of  the  formation 
between  40Kcal/mole  and  60Kcal/mole  are  within  the  transition  region. 


I-V  characteristic  of  sample  Nos.  3  and  8. 


5.5  Conclusions 


A  metal-amorphous  film-semiconductor  (MAS)  model  is  suggested 
to  be  the  physical  structural  model  of  a  real  Schottky  diode.  Based 
on  this  structural  model,  all  the  experimental  results  obtained  from 
GaP  (111)  Schottky  diodes  can  be  explained.  The  higher  Schottky 
barrier  heights  generally  obtained  by  C-V  methods  than  those  obtained 
by  other  methods  including  I-V  and  photo  emission  can  be  explained  by 
the  existence  of  the  amorphous  interlayer  in  the  diode  The  deviation 
of  diode  factor  from  its  ideal  value  can  be  attributed  to  this  amor¬ 
phous  film.  The  existence  of  the  amorphous  layer  has  been  demonstrated 
and  material  properties  such  as  the  resistivity  and  the  dielectric 
constant  of  the  amorphous  film  have  been  obtained  experimentally.  The 
structural  model  and  the  defect  model  are  combined  to  explain  the 
origin  of  Fermi-level  pinning  and  the  transition  of  interface  behavior 
between  the  covalent  and  ionic  semiconductors. 
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CHAPTER  6 

INVESTIGATION  OF  SURFACE- SCATTERING  LOSSES  OF  GaAs  WAVEGUIDE 

6.1  Introduction 

Integrated  optics  is  based  primarily  on  the  fact  that  light 
waves  can  propagate  through,  and  be  confined  by  very  thin  layers  of 
dielectric  materials.  All  optical  functions  must  be  obtained  in  a 
single  material  in  monolithic  integrated  optics.  These  optical 
functions  can  be  categorized  as  light  generation,  coupling,  modula¬ 
tion,  guiding  and  detection.  III-V  semiconductors  are  the  most 
promising  crystal  family  for  monolithic  integrated  optics  due  to 
their  high  electro-optic  coefficients  and  good  electrical  properties. 

These  materials  are  also  compatible  with  the  most  highly  advanced 
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fabrication  techniques,  including  epitaxial  growth  ,  diffusion  *  , 

84 

ion  implantation  and  ion-beam  sputtering. 

Two  approaches  can  be  used  to  achieve  the  integration  of  opti¬ 
cal  circuits.  One  is  to  fabricate  integrated  optical  components 
using  materials  best  suited  for  each  specific  device,  and  then 
integrating  all  of  the  components  into  a  circuit  using  an  appro¬ 
priate  substrate.  There  are  serious  problems  of  coupling  between 
various  components  using  this  technique,  and  it  is  difficult  to 
achieve  a  high  degree  of  integration.  The  other  approach  is  to 
fabricate  all  optical  components  on  a  single  substrate.  This 
approach  simplifies  the  problem  of  coupling  between  various 


127 


components  and  can  achieve  a  very  high  degree  of  integration.  There 
are,  however,  still  problems.  The  most  promising  substrate  material 
is  GaAs,  which  has  the  electrical  and  optical  properties  potentially 
suitable  for  achieving  high  level  integrated  optics.  The  optical 
transmission  losses  in  GaAs  still  cannot  be  reduced  to  an  acceptable 
level  (a  few  decibles  per  centimeter)  by  the  present  technologies  for 
material  preparation.  A  detailed  understanding  of  the  magnitude  of 
the  propagation  loss  in  semiconductor  waveguides  at  various  wave¬ 
lengths,  is  not  available.  It  becomes  important  to  gain  a  quanti¬ 
tative  understanding  of  the  surface  and  interface  properties  of  the 
III-V  compound  semiconductors  in  order  to  develop  a  fabrication 
technology  analogous  to  that  of  silicon  in  integrated  electronics. 

In  this  chapter,  the  results  of  a  quantitative  characterization 
of  the  surface  geometry  for  semiconductor  waveguides  is  reported. 

This  surface  geometric  information  will  be  correlated  to  formulate 
the  calculation  of  the  parameters  in  the  surface  scattering  theory 
and  to  obtain  theoretical  data  on  scattering  losses. 


6.2  Surface  Scattering  and  Losses 


The  foremost  barrier  in  integrated  optics  is  that  the  optical 
propagation  losses  in  III-V  compound  semiconductor  waveguides  cannot 
be  reduced  to  an  acceptable  level.  The  most  important  losses  in  a 
wave  propagation  along  a  thin  film  waveguide  are  (1)  scattering  from 
rough  surfaces  of  the  waveguide  during  the  Internal  reflection  along 
the  zigzag  path,  and  (2)  band-to-band  transitions.  Losses  caused  by 
the  first  mechanism  can  be  reduced  by  making  the  surface  of  the  wave¬ 
guide  as  smooth  as  possible.  Losses  caused  by  the  second  mechanism 
can  be  minimized  by  a  proper  choice  of  the  energy  gap  of  the  guiding 
material  relative  to  the  operating  wavelength. 
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As  reviewed  in  section  2.2.2,  a  theoretical  calculation  of  the 

surface  scattering  of  a  symmetric  slab  waveguide  was  reported  by 

Marcuse^  based  on  mode  conversions  to  the  guided  and  radiation 

modes.  His  results  needed  further  approximation  for  practical 

application.  A  simple  theory  of  the  surface  scattering  for  a  sym- 

18 

metric  waveguide  developed  by  Tien,  et  al.  was  introduced  in 

section  2.2.3  and  will  be  used  to  calculate  the  scattering  loss  of 

a  rib  waveguide.  Rib  waveguides  are  useful  for  lateral  confinement 
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and  inter-connections  in  monolithic  integrated  optics.  As  shown, 

15  -3 

Fig.  6.1  is  a  rib  waveguide,  comprised  of  a  lightly-doped  (10  cm  ) 
GaAs  film  of  height  t,  refractive  index  n^  with  a  rib  of  width  w  and 
height  Ah  etched  into  this  lightly-doped  film.  This  structure  is 
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sandwiched  between  a  highly-doped  (3x10  cm  )  GaAs  substrate  of 


refractive  index  and  air  of  refractive  index  n2>  where  both  n, 


Cross  section  (x,y)  and  top  view  (y,z)  of  a  rib 
waveguide,  the  effective  indices  N  and  Ng  of  the 
rib  and  the  surrounding  regions. 


Fig.  6.1 
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and  are  smaller  than  n^.  The  presence  of  free  carriers  in  a  semi¬ 
conductor  reduces  the  refractive  index  of  the  material.  This  effect 
is  sufficiently  strong  in  n-type  GaAs  to  enable  wave  confinement. 

The  depression  of  the  refractive  index  is  primarily  due  to  the  nega¬ 
tive  contribution  of  the  free-carrier  plasma  to  the  dielectric 
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constant.  This  contribution  for  n-type  GaAs  is  given  by 
An— 9.6xl0-21  N/nE2 - (6.1) 

where  N,n  and  E  are  the  free-carrier  concentration,  the  index  of 
refraction  of  pure  material,  and  the  energy  of  guiding  wave,  respec¬ 
tively.  The  rib  shape  provides  the  lateral  confinement,  and  will  be 
discussed  later. 

The  analysis  of  dielectric  rib  waveguides  is  much  more  complex 

chan  that  of  planar  guides,  and  no  exact  analytic  solutions  for  the 

modes  of  rib  waveguides  are  available.  A  simple  model,  the  effec- 
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tive  index  approximation  first  introduced  by  Knox  and  Toulios  ,  is 
followed  for  determining  the  characteristics  of  the  waveguide.  This 
method  is  closely  related  to  the  approximated  analytical  approach 
developed  by  Marcatili,  as  introduced  in  section  2.2.4. 

Local  solutions  for  the  y  dependence  are  solved  first  and  are 

used  to  determine  the  x  dependence  of  the  effective  refractive  index 

N  .  The  effective  index  model  attributes  the  effective  refractive 
e 

indices  to  the  rib  (r)  region  (Nf)  and  to  the  region  (s)  surrounding 

the  rib  (N  ) ,  as  3hown  in  the  top  view  of  the  Fig.  6.1. 
s 


The  effective  refractive  indices  describing  the  r  and  s  regions  are 
determined  by 


Ne-(8A/2) - (6.2) 

where  A  is  the  incident  wave  length  in  free  space.  The  propagation 
constant  8  in  the  uniform-layer  problem  is  defined  by 


8-  ✓  (nV-k2) - (6.3) 

where  k»2ir  \,n  is  the  refractive  index  in  the  guiding  layer,  and  k^ 
is  the  y  component  of  propagation  vector  in  the  guiding  layer. 

Taking  as  a  real  function  of  k^, 

S  *  /(n2-n b  k2-k2  - (6.4) 

i  i  y 

The  subscript  i  designates  the  layer  bounding  the  guiding  layer 
(Fig.  6.1).  i»l  denotes  the  substrate  layer  and  i«2  denotes  the 
cover  layer. 

There  are  two  families  of  modes  (E*  and  Ey  )  that  the  waveguide 

pq  pq  s 

can  support,  as  described  in  section  2.2.4.  For  the  E^j  mode,  the 
transdental  equation  is  similar  to  Eq.  (2.95), 

tan  k  t  »  k  (S.  +  S,)/(k2-S,  Sj - (6.5) 


This  can  be  used  to  evaluate  the  y  component  of  the  propagation 


vector  (ky)  within  a  guiding  layer  of  thickness  t.  The  effective 
refractive  indices  of  the  rib  region  (N^)  and  of  the  surrounding 
region  (Ng)  is  then  obtained  using  equations  (6.2)  and  (6.1). 

The  x  dependence  of  the  field  with  the  artificial  boundary  con¬ 
ditions  imposed  at  the  edges  of  the  rib  at  x*±l/2W  is  proportional  to 


cos  (kxX)  for  |X|<  1/2  W 

- (6.6) 

exp  (-|S3X|)  for  | X j >1/2  W 

where  S_,  a  real  function  of  k,  is  determined  by  0  ,0  and  k 

3»  i  r  s 

through  the  following  equation, 

S3  +  kx*0?’  “  8s - <6*7) 

where  k^  is  the  x  component  of  the  propagation  vector  and  can  be 
solved  using  transdental  equation  given  by  Eq.  (2.93);  e.g. 

tan  kW-  2  V3/((Ns/Nr)2  "  S3<VN3)2) - (6,8) 


A  computer  program,  (see  appendix) ,  can  be  used  to  obtain  kx  and  k 
Once  kx  and  k^  are  obtained,  the  propagation  constant  of  the  rib 
waveguide  can  be  obtained  using  Eq.  (2.91). 
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Consider  a  plane  wave  incident  on  the  upper  surface  of  the  rib 


covering  the  width  (W)  of  the  rib  in  the  y  direction,  as  shown  in 
Fig.  6.1.  To  cover  a  unit  length  of  the  rib  in  the  z  direction,  the 
plane  wave  has  a  width  of  cos0  in  the  direction  parallel  to  the  wave- 
front.  The  incident  power  P  carried  by  the  plane  wave  can  be  given 

P-(l/2)  nA2  W^eo/Uocose - (6.9) 

where  A  is  the  field  amplitude,  and  Eo  and  Uo  are  the  permittivity 
and  magnetic  permeability  of  free  space,  respectively.  The  power 
loss  Ap  at  two  interface  is  given  by 

Ap-(l/2)  nA2K2cos^0  W/£g /Ua  - (6.10) 

with  K-4tt/  (a2x  +  a22)1/Z 

where  a  .  and  a  ~  are  the  standard  deviation  of  the  interfaces 
rl  r2 

between  the  rib  and  cover  layers,  respectively.  The  total  power  P 
transmitted  in  the  z-direction  is  given  by 

P  »  n  sin9A2/4/y0/£o)  (t+Ah+l/s^+l/s^)  (w+ 

2(K2+s2)  /(s3(Ng/Nr)K2+s^(Nr/Ns)2)  - (6.11) 

The  power  attenuation  per  unit  length  of  the  rib  waveguide  by  means 
of  surface-scattering  losses  is  then  given  by 


ct  «  K2(cos29/sin9)  (t+Ah+l/s, +1/33)  *  (w+ 
2(K2+s2)  /  <s3 (Ng/Nr)  V+s2  (Nr/Ng) V1 


(6.12) 


The  first  term  K  in  the  equation  depends  only  on  the  surface  proper¬ 
ties  of  interfaces.  The  experimental  techniques  used  for  investiga¬ 
ting  the  surface  properties  of  interfaces  will  be  discussed  in  the 


next  section. 
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6.3  Outline  of  Experiments 

Rib  waveguide,  described  in  the  previous  section,  was  used  to 
investigate  the  surface-scattering  losses  of  III-V  compound  semi¬ 
conductor  waveguides.  The  waveguide  can  be  prepared  on  highly-doped 
n-type  GaAs  crystals.  Both  (111)  and  (100)  orientations  of  GaAs 
wafers  were  used.  The  doping  concentrations  of  these  wafers  were 

determined  by  resistivity  measurements  which  were  carried  out  by  the 
27  28  29 

four-point  probe  ’  and  Van-der  Pauw  methods,  as  stated  in  sec¬ 
tion  3.4.1.  The  GaAs  crystals  are  1.5  in. -diameter  by  0.02  in. -thick 
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wafers,  and  chemo-mechanically  polished  with  Q.075Z  bromine- in¬ 
methanol  solution  using  a  PAN-W  polishing  cloth.  Five  samples  were 

cut  from  the  same  wafer  by  cleaving.  These  samples  with  surface 

2 

areas  varying  from  about  2.S  to  3  cm  were  cleaned  in  an  ultrasonic 
cleaner  using  trichloroethylene,  acetone,  methonal  and  D. I. -water, 
and  then  attached  to  a  sample  holder  using  Apiezon  wax.  Different 
degrees  of  surface  roughness  were  obtained  by  polishing  the  samples 
with  alumina  powders  of  different  particle  sizes.  Five  different 
sizes,  0.05,  0.3,  0.5,  1.0  and  3.0Vms,  were  used.  A  new  polishing 
cloth  was  Installed  once  the  size  of  polishing  powders  was  changed. 

The  speed  of  the  polisher  was  set  at  a  fixed  rate  for  each  polishing 
process.  The  sample  holder  was  manually  controlled  to  provide  a 
random  motion.  Distilled  water  was  used  in  the  polishing  process. 

After  4  minutes  of  polishing,  the  samples  were  removed  from  the 
sample  holder  using  trichloroethylene  and  rinsed  in  an  ultrasonic 
cleaner  using  acetone,  methanol  and  D. I. -water.  The  samples  were 
blown  dry  with  nitrogen  and  ready  for  surface  roughness  analysis. 


Laser  scattering  techniques  ,  described  in  section  3.2,  were 

performed  to  characterize  the  surface  geometry  of  the  GaAs  substrates. 

The  incident  laser  beam  of  0.5mm  beam  diameter  scanned  over  an  area  of 
2 

9mm  around  the  center  of  each  sample.  Five  locations  of  each  sample 
shown  in  Fig.  3.4  were  first  rotationally  scanned  by  the  incident 
beam.  The  directions  at  which  the  maximum  and  minimum  scatter  occurred 
were  determined.  Each  sample  was  then  scanned  translationally  along 
these  two  directions  separately.  The  average  BRDF  value  by  averaging 
the  BRDF  values  in  these  two  directions  was  obtained.  The  average  BRDF 
value  quantitatively  characterizes  the  degree  of  surface  roughness. 

The  last  stage  of  Che  work  is  the  fabrication  of  the  rib  wave¬ 
guide  and  the  study  of  the  optical  propagation  losses  of  such  wave¬ 
guides.  The  first  step  of  this  stage  is  to  grow  a  lightly-doped 
(~10A  cm  )  GaAs  thin  film  (~2y)  upon  the  substrates.  In  order  to 
preserve  the  original  surface  morphology  of  the  substrates,  the  only 
possible  growth  technique  can  be  used  is  molecular  beam  epitaxy  (MBE) . 
The  growth  race  of  this  technique  can  be  controlled  in  the  molecular 
order  such  that  the  surface  morphology  of  the  substrates  can  be  kept. 

The  next  step  is  to  form  the  rib  on  the  lightly-doped  thin  film.  As 

25  90  91 

described  In  section  3.3,  anodic  oxidation  ’  *  can  be  used  in 

forming  the  rib  geometry  with  a  rib  of  width  (w«3y)  and  height 
(Ah-.2y)  etch  into  the  thin  film.  Then  the  optical  propagation 
losses  of  such  waveguides  will  be  measured  and  be  used  to  develop  the 
correlation  between  the  surface  roughness  obtained  from  laser  scat¬ 
tering  measurements  and  surface-scattering  parameters  or  ^ >ar2  an<*  K* 
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6.4  Results  and  Discussion 

The  carrier  concentrations  of  GaAs  wafers  in  (100)  and  (111) 

18 

directions  determined  by  resistivity  measurements  were  2.4+  .2x10 
—3  18  —3 

cm  and  2.1+0.3x10  cm  ,  respectively.  According  to  equation  6.1, 

the  depressions  of  the  refractive  index  from  that  of  the  pure  material 

-3  -3 

at  E»1.24ev  were  -4.5x10  and  -3.7x10  for  (100)  and  (111)  orienta¬ 
tions,  respectively.  The  refractive  index  change  between  pure  and 
the  n-type  wafers  are  large  enough  to  confine  the  light  vertically  if 
a  pure  or  lightly-doped  layer  of  sufficient  thickness  can  be  fabri¬ 
cated  on  the  n-type  GaAs  substrates. 

10  samples  were  quantitatively  analyzed  by  laser  scattering 
technique.  They  were  labeled  in  accordance  to  the  crystal  orientation 
of  the  substrate  and  the  grit  size  of  the  polishing  powders  which  were 
used  to  generate  the  surface  roughness.  As  stated  in  Section  3.2, 
five  locations  of  each  sample  depicted  in  Fig.  3.4  were  first  rota- 
tionally  scanned  by  the  incident  beam.  A  circular  plot  of  scatter 
intensity  vs.  rotation  angle  was  obtained  for  each  sample,  as  shown 
in  Fig.  6.2.  The  curve  provided  a  visual  impression  of  the  directions 
at  which  the  maximum  and  minimum  scatter  occurred.  Each  sample  was 
then  translationally  scanned  by  the  incident  beam  along  the  maximum 
and  minimum  scattering  directions  separately.  In  Fig.  6.3,  a  three- 
dimensional  plot  shows  the  BRDF  values  of  the  sample  scanned  by  the 
incident  beam  along  one  of  the  two  directions.  As  shown  in  Fig.  6.4, 
average  plots  were  made  by  averaging  perpendicular  to  the  direction 
of  scan,  and  an  average  line  was  obtained.  Finally,  the  overall 


i 


Average  BRDF  values  obtained  by  averaging  perpendicular 
to  the  direction  of  scan  where  (a)  is  along  the  minimum 
scatter  direction,  (b)  is  along  the  maximum  direction 
and  (c)  is  the  average  of  (a)  and  (b) . 
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average  scattered  number  (BRDF)  as  listed  in  Table  6.1  was  obtained 
which  quantitatively  represented  the  degree  of  surface  roughness  of 
each  sample.  Samples  1  through  5  and  samples  6  through  10  were 
cleaved  from  the  same  GaAs(lOO)  and  GaAs(lll)  wafers,  respectively. 
They  were  polished  through  the  same  procedure.  The  only  different 
parameter  in  the  entire  procedure  between  samples  1  through  5  was 
the  size  of  powders  used  to  polish  them.  The  overall  average  scat¬ 
tered  number  (average  BRDF)  experimentally  obtained  for  each  sample 
was  plotted  against  the  size  of  the  powders  used  to  polish  it  as 
shown  in  Fig.  6.5.  Similar  plot  for  samples  6  through  10  was  shown 
in  Fig.  6.6.  Both  figures  showed  that  the  degrees  of  surface  rough¬ 
ness  of  the  samples  were  proportional  to  the  grit  size  of  the  powder. 
The  linear  correlation  between  grit  sizes  shows  that  laser  scattering 
technique  can  give  a  quantitative  characterization  of  the  degrees  of 
surface  roughness. 

The  last  stage  of  the  experiment  is  currently  on  going.  Further 
work  is  needed  in  order  to  complete  the  measurements  of  the  signal 
losses  of  the  waveguides  with  different  degrees  of  surface  roughness, 
the  development  of  the  correlation  between  the  surface  roughness  and 
the  surface-scattering  parameters  and  finally,  the  investigation  of 
the  optical  properties  of  the  waveguides. 


Table  6.1 


Surface  Roughness  of  Samples  1  Through  10 


Sample  No.  Label  Orientation  Grit  Size(  m)  Ave  BRDF(PPM/ST) 


1 

GAl 00005 

100 

.05 

40.61 

2 

GA10003 

100 

.3 

2532.28 

3 

GAl 0005 

100 

.5 

9235.28 

4 

GAl 00 10 

100 

1.0 

49523.00 

5 

GAl 0030 

100 

3.0 

352.89 

6 

GA11103 

111 

.3 

1170.61 

7 

GAl 1105 

111 

.5 

12114.82 

8 

GAl 1110A 

111 

1.0 

15866  00 

9 

GA11110B 

111 

1.0 

41219.00 

10 

GAl 11 30 

111 

3.0 

3205.56 

SURFACE  ROUGHNESS  (PPM/STR  ) 

Fig.  6.5  The  overall  average  number  (surface  roughness)  plotted  vs.  the 
grit  size  of  polishing  powders  for  samples  No.  1  through  No.  5 


GaAs  (til) 


SURFACE  ROUGHNESS  (pPM/STR) 

The  overall  average  number  of  samples  No.  6  through  No.  10 
plotted  vs.  the  grit  size  (pm)  of  polishing  powders. 


CHAPTER  7 


SUMMARY  AND  CONCLUSIONS 

The  cheorecical  formulation  for  determination  of  the  atomic 
structure  of  solid  surface  by  LEED  intensity  analysis  was  reviewed. 

A  dynamical  multiple  scattering  computer  program  provided  by 
S.  Y.  Tong,  et  al.  was  slightly  modified  for  calculating  the  theo¬ 
retical  I-V  curves  of  GaP  (110)  face.  The  intralayer  multiple 
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scattering  was  treated  exactly  using  direct  summation  method  while 
for  the  interlayer  multiple  scattering,  the  renormalized  -  forward  - 
scattering  (RFS)^»^®  perturbation  theory  was  used.  The  inner  poten¬ 
tial  was  taken  to  be  8eV  and  the  inelastic  damping  was  set  equal  to 
SeV.  The  electron-phonon  interaction  has  been  incorporated  into  the 
theoretical  formalism.  Both  room  and  low  temperatures  LEED  experi¬ 
mental  I-V  curves  were  recorded  and  essentially  identical.  I-V 
curves  obtained  for  the  first  16  inequivalent  nonspecular  beams 
emerging  from  the  sample  were  compared  with  the  theoretical  curves 
calculated  from  17  different  proposed  structures.  Our  study  has 
indicated  that  the  atomic  structure  of  GaP  (110)  face  is  quite  simi¬ 
lar  to  that  of  GaAs  (110)  regarding  the  ripple  geometry.  The  P 

atoms  are  rotated  outward  from  the  GaP  (110)  plane  and  Ga  atoms  are 

o 

rotated  inward  from  the  plane  by  an  angle  of  27  in  order  to  preserve 
the  bond  lengths  between  the  nearest-neighbor.  Furthermore,  the  top 
most  layer  is  compressed  by  5Z  such  that  the  top  layer  spacing  is 

9 

reduced  by  approximately  0.1A.  The  experimental  LEED  patterns  of 
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GaP  (111)  face  were  studied.  It  is  concluded  that  a  good  Schottky 


contact  should  be  relatively  easy  to  obtain  on  the  GaP  (111)  face. 

The  complete  surface  structure  analysis  of  GaP  (111)  is  premature 
and  needs  further  study. 

As  suggested  in  the  previous  LEED  studies,  Schottky  barriers 
have  been  fabricated  on  n-type  GaP  (111)  crystals  to  study  the  basic 
physical  mechanism  for  the  metal-semiconductor  interaction  and  for¬ 
mation  of  Schottky  barriers.  Capacitance-Voltage  (C-V) ,  capacitance- 
frequency  (C-f) ,  current-voltage  (I-V)  and  conductance-frequency 
were  performed  to  analyze  the  electrical  properties  of  the  Schottky 
diodes.  AES  studies  were  used  to  study  the  interface  chemical  com¬ 
positions.  Strong  interdiffusion  of  the  metal  and  semiconductor 
atoms  at  the  interface  region  was  observed  in  the  AES  studies  of 
Au-GaP  and  Al-GaP  Schottky  barriers.  Based  on  the  experimental 

observations,  a  metal-amorphous  film-semiconductor  (MAS)  configura- 
72 

tion  is  suggested  to  be  the  physical  structural  model  of  a  real 

Schottky  barrier.  The  composition  of  the  amorphous  film  is  a  random 

mixture  of  metal  and  semiconductor  atoms.  The  thickness  of  this 

film,  W,  determined  experimentally  from  AES  depth-composition  pro- 

0  0 

files,  is  in  the  range  of  50A  to  500A  depending  upon  fabrication 
details.  An  electrical  equivalent  model  is  derived  on  the  basis  of 
the  structural  model  to  describe  the  fundamental  behaviors  of  a  real 
Schottky  diode.  The  deviation  of  the  diode  factor  n  from  its  ideal 
value  of  unity  is  usually  observed,  and  it  has  been  demonstrated 
theoretically  that  a  higher  value  of  diode  factor  can  be  attributed 
to  the  existence  of  the  amorphous  film.  The  material  properties 
such  as  the  resistivity  and  the  dielectric  constant  of  the  amorphous 


film  have  been  obtained  experimentally.  The  model  has  provided  a 
reasonable  explanation  for  the  higher  Schottky  barrier  heights 
generally  reported  by  C-V  measurement  than  those  obtained  by  other 
methods  including  I-V  and  photo  emission.  Based  on  the  MAS  struc¬ 
tural  model,  extrinsic  surface  states  on  the  semiconductor  would  be 

created  before  an  amorphous  film  is  formed,  which  strongly  supports 
68  69 

the  defect  model  *  proposed  to  explain  the  Fermi-level  pinning  of 
a  real  Schottky  barrier.  Our  MAS  model  can  also  be  combined  with 
the  defect  model  to  explain  the  transition  of  interface  electronic 
behavior  between  covalent  and  ionic  semiconductors.  However,  our 
AES  result  cannot  rule  out  the  probability  of  Fermi-level  pinning 
caused  by  a  contamination  of  1%  of  a  monolayer  at  the  metal- 
semiconductor  interface  due  to  the  limitation  of  the  Auger  sensi¬ 
tivity.  Further  study,  in  an  interlock  system,  which  has  the 
capability  of  fabricating  Schottky  diodes  in  an  UHV  chamber  on  a 
newly  cleaved  or  ion-sputtering  cleaned  semiconductor  surfaces  to 
eliminate  such  contamination,  is  needed  and  currently  underdeveloped. 
Studies  involving  different  metals  and  other  III-V  compound  semi¬ 
conductors  may  be  used  in  the  future  to  demonstrate  the  suggested 
physical  structural  model  of  a  real  Schottky  barrier. 

The  theoretical  calculation  of  surface  scattering  theory  was 

18 

reviewed.  The  simple  theory  of  surface  scattering  and  the  effec- 

87 

tive  index  approximation  were  used  to  calculate  the  surface¬ 
scattering  losses  of  a  rib  waveguide.  The  doping  concentration  of 

the  GaAs  substrates  in  (100)  and  (111)  directions  obtained  by 

18  -3  18  —3 

resistivity  measurements  were  2.4x10  cm  and  2.1x10  cm  , 
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respectively.  The  vertical  confinement  of  the  light  in  the  waveguide 
can  be  provided  if  a  lightly-doped  layer  of  sufficient  thickness  is 
fabricated  on  the  highly-doped  GaAs  substrates.  Different  degrees 
of  surface  roughness  were  prepared  by  polishing  the  samples  with 
alumina  powders  of  different  particle  sizes.  Laser  scattering 
technique  was  then  conducted  to  characterize  the  surface  geometry  of 
these  samples.  The  overall  average  scattered  numbers  (average  BRDF) 
obtained  for  all  the  samples  showed  that  the  degrees  of  surface 
roughness  of  the  samples  were  proportional  to  the  grit  size  used  to 
polish  them.  We  conclude  that  the  laser  scattering  method  can  be 
used  quantitatively  to  characterize  the  degrees  of  surface  roughness. 

Further  work,  comparing  the  experimental  measurements  of  the 
signal  losses  obtained  from  the  waveguides  with  different  degrees 

of  surface  smoothness  with  the  theoretical  data  of  scattering  loss, 

/ 

is  needed  in  order  to  develop  the  correlation  between  the  surface 
roughness  and  the  parameters  in  the  surface  scattering  theory;  and 
finally,  to  investigate  the  optical  properties  of  the  waveguides. 
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APPENDIX 


Calculation  of  the  Propagation  Constant  in  a 
Rectangular  Waveguide 

The  following  programs  calculate  the  propagation  constant  6  of 

a  waveguide  propagating  along  z-direction  in  a  rectangular  dielectric 

waveguide  with  dimension  a  and  b  as  shown  in  Fig.  2.9.  The  program 

uses  Newton's  method  to  solve  the  transdental  equations  (2.93)  and 

(2.95),  for  k  and  k  ,  respectively.  Once  k  and  k  are  determined, 
x  y  x  y 

the  propagation  constant  8  can  be  obtained  from  Eq.  (2.97).  Program¬ 
ming  was  in  standard  fortran  IV,  and  programs  have  been  tested  in  an 
IBM  360  under  WATFTV  compilers. 

The  first  step  calculates  k  and  k  .  Eqs.  (2.93)  and  (2.95) 

x  y 

can  be  written  as 

2  2  2  2 

f  (kx)*kxa-pu— arctanCn^s^/Ojk^-arctan^s^/n^k^-—  (A-I) 


and 


f(k  )»k  b-qir-arctan(s-/k  )-arctan(s. /k  ) - - (A-2) 

yy  Zy  4  y 

where  the  integers  p  and  q  indicate  the  number  of  maxima  of  the  field 
distribution  in  x  and  y  directions.  Newton's  method  for  finding  zeros 
of  f(x)  is  most  widely  known,  and  it  can  be  expressed  by  an  iteration 


process 


\+i  ■  vc(xv)/r(\) 


(A-3) 
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Afcer  k  and  k  are  obcained, 
x  y 


che  propagaclon  constant  3  can  be 


determined . 


1 


COMMENT  PROGRAM  TO  CALCULATE  THE  PROPAGATION  CONSTANT  IN  A 
C  RECTANGULAR  WAVEGUIDE  USING  NEWTON'S  METHOD 

COMPLEX  CN1 , CXX , CS  3 , CS  5 , CDS  3 , CDS  5 , CFKX , CFKY 
COMPLEX  CSQRT , CIN , CEN , CS2 , CS4 , CDS 2 ,CDS4 , CN2 
COMPLEX  CN3 , CN4 , CN5 , CBETA , CMPLX , CATAN 
PHI-3.14159625 
CN1 -CMPLX ( 3 . 4700 , -0 . 0005 ) 

CN2-1.0 

CN3-3.42 

CN4-3.4 

CN5-3.42 

RK-2*PHI*1. 424/1. 44 

A-3.0 

B-1.0 

C  THIS  SECTION  CALCULATES  CKX 

N-l 
P-1.0 

CKX-CMPLX(0. 4027, -0.005) 

15  CS3-CSQRT((CN1**2-CN3**2)*RK**2-CKX**2) 

CS5-CSQRT ( (CN1**2-CN5**2) *RK**2-CKX**2) 

CDS3- (-1 ) *CKX*CN3**2/CN1 **2/CS3 
CDS5-(-l)*CKX*CN5**2/CNl**2/CS5 
CFKX«CXX*A-P*PHI+CATAN ( -CDS  3 )  +CATAN ( -CDS  5 ) 
CDFKX-A+(CS3*CN**2/CN1**2-CDS3*CKX) / ( 1+CDS3**2) 
1/CS3**2+(CS5*CN5**2/CN1**2-CDS5*CKX) / (1+CDS5**2) 


1/CS5**2 


CIN-CFKX/CDFKX 
CKXN-CKX-CIN 
WRITE ( 6 , 4)  N , CKXN 

4  FORMAT (4H,  N-  ,14, 7H,  CKXN-  ,2E19.9) 

IF(CABS(CIN)-.1E-07)1Q,17 ,17 
17  CKX-CKXN 

IF (N-200) 21,22,22 

21  N-NH 

GO  TO  15 

22  WRITE (6 ,23) CKX 

23  FORMAT (16H,  CHOOSE  NEW  CKX.2E19.9) 

10  WRITE (6, 24) N, CKX 

24  FORMAT (4H,  N-  ,14,6H,  CKX-  ,2E19.9) 

COMMENT  THIS  SECTIONCALCULATE  CKT 

M-l 
Q-l  .0 

CKY-( 1.8334, -0.002) 

45  CS2-CSQRT( (CN1**2-CN2**2)*RK**2-CKY**2) 

CS4-CSQRT ( (CN1**2-CN4**2) *RK**2-CKT**2) 
CDS2-(-l)*CKY/CS2 
CDS4-(-i)*CKY/CS4 

CFKY-CKY*B-Q*PHI+CATAN ( -CDS  2 ) +CATAN ( -CDS4 ) 
CDFKY+B+<CS2-CDS2*CKT) / (1+CDS2**2) /CS2**2 
1+(CS4-CDS4*CKY) / ( 1+CDS4**2) /CS4**2 
CEN-CFKY / CDFKY 
CKYN-CKY-CEN 
WRITE(6,44)M,CKTN 


44 


FORMAT (4H,  M-  •,14,7H,  CKYN-  .2E19.9) 
IF(CABS(CEN)-.lE-07)40,40,47 
47  CKY-CKYN 

IF(M-200)41 ,42,42 

41  M-Mfl 

GO  TO  45 

42  WRITE(6,43)CKY 

43  FORMAT (16H,  CHOOSE  NEW  CKY  .2E19.9) 

40  WRITE (6 ,64) M, CRY 

64  FORMAT (4H,  M-  ,14,6H,  CKY-  .2E19.9) 

COMMENT  THIS  SECTION  CALCULATE  CBETA 

CBETA-CSQRT ( CN 1 **2*RR**2-CKX**2--CKY** 2 ) 

WRITE (6, 65) CBETA 

65  FORMAT (8H,  CBETA-  ,2E19.9) 

STOP 

END 

C 

COMMENT  FUNCTION  CATAN  PROVIDES  THE  VALUE  OF  COMPLEX  ARCTANGENT 
FUNCTION  CATAN (CZ) 

COMPLEX  CATAN, CZ.CMPLX, Cl 
CI-CMPLX(O.O.O.l) 

X-REAL(CZ) 

Y-AIMAG(CZ) 

CATAN-ATAN ( 2*X/ ( 1 -X**  2-Y**2) ) / 2+CI /4*AL0G( ( 1+Y) **2 
1+X**2) / ( (l-y)**2+X**2) ) 


RETURN 


